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1  INTRODUCTION 


Computational  electromagnetics  continues  to  rely  heavily  on  surface  integral  equations  for 
the  efficient  numerical  solution  to  scattering  from  perfectly  electrically  conducting  (PEC) 
bodies  [1],  [2],  [3].  Both  the  magnetic-field  integral  equation  (MFIE),  which  was  derived  by 
Murray  [4]  in  1931,  and  the  electric-field  integral  equation  (EFIE),  which  was  derived  along 
with  the  MFIE  in  the  definitive  1949  paper  by  Maue  [5],  are  applied  only  to  the  surface 
current  of  the  scatterer  and  thus  require  a  number  of  unknowns  proportional  to  the  sur¬ 
face  area  in  square  wavelengths  of  a  three-dimensional  (3-D)  scatterer.^  The  MFIE  is  an 
absolutely  convergent  integral  equation  of  the  second  kind  whose  associated  discrete  matrix 
has  a  bounded  condition  number  as  the  number  of  discretizations  (grid  points)  on  the  sur¬ 
face  of  the  scatterer  approaches  an  infinite  value.  The  EFIE  is  a  conditionally  convergent 
integro-differential  equation  of  the  first  kind  whose  associated  discrete  matrix  has  an  un¬ 
bounded  condition  number  as  the  number  of  discretizations  on  the  surface  of  the  scatterer 
approaches  an  infinite  value  [5],  [6],  [7],  [8].  Therefore,  the  determination  of  surface  current 
(and  subsequently  the  scattered  fields)  for  a  scatterer  using  the  EFIE  generally  requires  a 
more  sophisticated  numerical  scheme  with  more  unknowns  per  square  wavelength  than  the 
MFIE.  However,  the  matrix  solution  to  the  MFIE  degenerates  to  an  underdetermined  set 
of  equations  for  the  surface  current  on  an  open  scatterer  (infinitesimally  thin  conductor)  [9, 
pp.  168-172]  and  approaches  degeneracy  if  parts  of  the  conducting  surface  become  much 

closer  to  each  other  than  a  wavelength  as  in  the  case  of  a  thin  plate,  thin  wire,  or  tip  of  a 
^Although  we  often  refer  to  “the”  MFIE  or  “the”  EFIE,  these  equations  can  be  written  in  many  different 

forms. 
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small-angle  cone.  Consequently,  to  determine  the  scattering  from  many  different  kinds  of 
PEC  bodies  using  surface  integral  equations,  the  EFIE  must  be  employed. 

Both  the  EFIE  and  the  MFIE  have  a  serious  limitation.  As  Murray  [4]  and  Maue  [5] 
pointed  out,  the  EFIE  and  the  MFIE  fail  to  produce  a  unique  solution  for  the  current  on  a 
PEC  scatterer  at  frequencies  equal  to  the  resonant  frequencies  of  the  interior  cavity  formed 
by  the  surface  of  the  scatterer.  Since  the  density  of  the  cavity  resonant  frequencies  increases 
rapidly  with  frequency  beyond  the  first  resonance,  which  for  many  cavities  occurs  when  the 
maximum  linear  dimension  is  approximately  one  wavelength,  the  numerical  solution  of  3-D 
multiwavelength  bodies  is  severely  hampered  by  these  spurious  resonances.  In  principle,  the 
spurious  resonances  in  the  solution  for  the  surface  current  occur  only  at  discrete  frequencies, 
and  for  the  EFIE  they  should  not  contribute  to  the  scattered  field  [6].  In  numerical  practice, 
however,  they  contaminate  the  surface  current  and  scattered  fields  of  both  the  EFIE  and 
MFIE  over  finite  bandwidths  about  the  cavity  resonant  frequencies.  Within  these  band- 
widths  the  determinants  of  the  solution  matrices  of  the  integral  equations  become  too  small, 
and  the  condition  numbers  become  too  large,  to  produce  an  accurate  numerical  solution. 

Among  the  several  ways  to  eliminate  the  spurious  resonances  from  the  original  integral 
equations,  the  combined-field  integral  equation  (CFIE)  and  the  dual-surface  magnetic-field 
integral  equations  (DSMFIE)  have  distinguished  themselves  in  numerical  practice  as  ef¬ 
fective,  efficient,  and  convenient  alternative  integral  equations  for  eliminating  the  spurious 
resonances  [9,  pp.  222-225],  [10],  [11],  [12],  [13],  [7],  [14],  [2,  ch.  6]. 

The  CFIE  [15],  [9],  [10],  [11]  is  formed  by  combining  the  MFIE 

n  X  H*""(r)  =  J(r)/2  -nx  J  J(r')  x  V'G(r,  T')dS'  (1) 

s 


2 


with  the  EFIE 


to  get 


n  X  E*”"(r)  =  -^n  x  f  [k'^3{r')G{r,T')  -  V'5  •  J(r')V'G'(r,  r')]rf5' 

(jJCci  J 


n  X  (1  -  ao)H*”=(r)  -  x  F""(r) 

Zo  J 

=  (1  -  ao)  J{r)/2  -  n  X  I  J(rO  x  VG{T,r')dS' 

s 

-^n  X  n  X  y [A:^ J(r')G'(r,  r')  -  V'5  •  J(r')  V'C?(r,  T')]dS'  (3) 

s 

where  J(r)  is  the  surface  current  density  on  the  surface  S  of  the  scatterer  (r  €  5),  Vs*  is 
the  surface  divergence  [16,  Appendix  2,  18.],  n  is  the  outward  unit  normal  from  the  surface 
5,  eo,  /iO)  and  Zq  =  y^/co  are  the  permittivity,  permeability,  and  impedance  of  free  space, 
and  are  the  electric  and  magnetic  fields  incident  upon  the  scatterer.  Harmonic 

time  dependence  of  the  form  exp{j(jjt)  has  been  suppressed,  k  =  and  the  free-space 

Green’s  function  is  given  by 


exp(-iA:|r  -  r'|) 
4xlr-r'l  ■ 


The  combination  parameter  ao  is  a  real  constant  that  can  assume  values  between  0  and  1, 
with  ao  =  .2  as  a  typical  choice  [17].  Mautz  and  Harrington  [10],  [11]  give  a  straightforward 
proof  that  the  CFIE  has  a  unique  solution  at  all  frequencies  and  therefore  eliminates  the 
spurious  resonances  from  the  original  MFIE  and  EFIE. 

The  DSMFIE  and  dual  surface  electric-field  integral  equation  (DSEFIE)  can  be  written 
as  [13],  [7] 

n  X  H‘^(r)  =  J(r)/2  -  n  x  |  J(r')  x  VG^{r,  r')d5"  (5) 

s 


3 


(6) 


n  X  E‘'(r)  =  -^n  x  [ [A:2j(r')C?‘'(r,r')  -  Vg  •  J(r')V'C?‘^(r,  r')]d5' 

(jJCq  y 

where  H‘^(r),  E‘^(r),  and  G‘^(r,  r')  are  defined  as 

H‘'(r)  =  H^"=(r)  +  (7) 

Ed(r)  =  E*""(r)  +  aE*”"(r‘')  (8) 

G‘'(r,r')  =  G(r,r')  +  QG'(r^r')  (9) 

with  r  G  S'  and  lying  on  a  surface  S*^  a  normal  distance  S  inside  S  as  shown  in  Fig.  1. 

The  DSMFIE  and  DSEFIE  in  (5)  and  (6),  although  identical  in  form  and  comparable  in 
complexity  to  the  original  MFIE  and  EFIE,  each  provide  a  unique  solution  for  J  (r)  at  all  real 
frequencies  as  long  as  the  combination  constant  a  has  an  imaginary  part  and  the  positive 
real  separation  constant  <5  is  less  than  half  a  wavelength  (A/2)  [13].  In  recent  numerical 
practice,  we  usually  choose  a  equal  to  .25/  and  5  between  A/32  and  A/4. 

It  is  not  necessary  for  uniqueness  of  solution  to  keep  the  inner  surface  at  a  fixed  distance 
5  from  the  scattering  surface  as  long  as  ^  is  less  than  A/2  and  large  enough  to  numerically 
maintain  uniqueness  of  solution.  In  practice  it  is  not  always  necessary  to  include  the  inner 
surface  for  every  observation  point  or  to  make  the  inner  surface  a  continuous  surface.  This 
means  that  one  can  often  avoid  difficulties  in  defining  the  inner  surface  for  scatterers  with 
complicated  shapes  by  simply  leaving  out  the  inner  surface  if  it  gets  closer  than  a  minimum 
distance  5  from  any  part  of  the  surface  S. 

The  dual-surface  integral  equations  were  introduced  a  number  of  years  ago  [18]  and  the 
DSMFIE  (5)  has  been  successfully  applied  to  multiwavelength  rectangular  boxes  [13],  [7]  and 
to  bodies  of  revolution  [14],  [19].  To  date,  a  numerical  solution  to  the  DSEFIE  (6)  for  3-D 
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Figure  1.  Geometry  of  a  perfectly  conducting  scatterer  with  surface  S  and  a  fictitious  “dual 
surface”  5*^  at  a  normal  distance  S  inside  S.  (The  value  of  6  can  vary  over  S.) 
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scatterers  has  not  been  successfully  formulated  and  programmed^,  and  thus  it  is  the  primary 
aim  of  this  report  to  formulate,  test,  and  document  the  DSEFIE  (6)  for  bodies  of  revolution 
(BOR’s). 

What  is  the  utility  of  a  dual-surface  integral  equation  computer  program  for  BOR’s 
when  CFIE  computer  programs  for  BOR’s  are  available  [21]?  We  can  answer  this  ques¬ 
tion  by  reviewing  some  of  our  recent  work  with  surface  integral  equations  in  connection 
with  incremental  length  diffraction  coefficients  (ILDC’s).  In  order  to  test  the  accuracy  of 
shadow-boundary  ILDC’s  developed  to  improve  upon  the  accuracy  of  physical  optics  (PO) 
for  calculating  high-frequency  scattering,  we  applied  them  to  multiwavelength  perfectly  con¬ 
ducting  BOR’s  and,  in  particular,  to  prolate  spheroids  [19].  To  determine  the  accuracy  of 
the  calculated  PO-fILDC  solution,  we  needed  a  highly  accurate  (<  ±.l  dB  error  in  the 
far-field  patterns)  numerical  solution  to  the  exact  field  equations  for  scattering  from  prolate 
spheroids.  (No  eigenfunction  solution  exists  for  electromagnetic  scattering  from  spheroids.) 
We  used  our  in-house  DSMFIE  BOR  program  to  compute  the  scattering  from  the  prolate 
spheroids.  (As  explained  above,  the  dual  surface  is  required  to  eliminate  the  spurious  reso¬ 
nances  from  the  solution  to  the  MFIE.)  Although  we  gained  some  confidence  in  the  DSMFIE 
solution  by  noting  convergence  of  the  solution  as  the  surfaces  of  the  spheroids  were  divided 
into  smaller  and  smaller  increments,  and  by  comparing  the  computed  and  exact  Mie  series 

solutions  for  the  sphere,  we  could  only  be  sure  that  the  DSMFIE  numerical  solution  was  con- 
method  of  moments  solution  for  the  DSEFIE  applied  to  BOR’s  was  considered  in  [20].  Although 

the  details  of  the  method  of  moments  formulation  are  not  given  in  [20],  they  can  be  determined  from  the 

computer  program  written  in  conjunction  with  [20].  An  analysis  of  this  method  of  moments  formulation 

revealed  that  it  was  invalid  for  the  dual  surface  part  of  the  DSEFIE. 
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verging  accurately  enough  to  the  correct  solution  for  eccentric  spheroids  if  we  could  compare 
it  to  another  independent  numerical  solution. 

To  obtain  a  second  independent  numerical  solution  for  scattering  from  the  prolate  spheroids, 
we  used  a  CFIE  computer  program  called  CICERO  [21] .  The  CFIE  and  the  DSMFIE  agreed 
over  the  entire  far-field  pattern  to  within  the  thickness  of  the  lines  used  to  plot  the  fax-field 
patterns.  Even  though  the  CFIE  and  DSMFIE  are  each  integral  equations,  the  CFIE  com¬ 
bines  both  the  electric  and  magnetic  field  operators  to  form  a  surface  integral  equation  that 
is  very  different  from  the  DSMFIE,  which  involves  only  the  operator  of  the  magnetic  field 
integral  equation.  Therefore,  the  close  agreement  of  the  two  very  different  integral  equation 
solutions  gave  us  the  confidence  to  designate  these  solutions  as  the  “highly  accurate  reference 
solution”  to  which  we  could  compare  our  approximate  high-frequency  ILDC  solution. 

This  experience  not  only  confirmed  in  us  an  appreciation  for  surface  integral  equations, 
but  it  also  clearly  demonstrated  the  importance  of  having  available  two  independent  sur¬ 
face  integral  equation  formulations  that  can  be  applied  and  solved  straightforwardly  and 
efficiently.  In  the  work  with  spheroids  described  above,  the  two  surface  integral  equations 
we  used  were  the  DSMFIE  and  the  CFIE.  However,  as  explained  above,  for  infinitesimally 
thin  conductors  (open  surfaces),  the  magnetic  field  operator  degenerates  to  give  an  under¬ 
determined  set  of  equations,  and  for  conductors  containing  narrow-angle  wedges  and  tips, 
the  magnetic  field  operator  becomes  unstable.  Thus,  for  many  scattering  geometries,  such 
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as  the  small-angle  cone-sphere,  one  must  use  the  DSEFIE  instead  of  the  DSMFIE.^ 

Consider  the  following  example  of  using  the  DSEFIE  where  the  DSMFIE  has  difficulty 
producing  as  accurate  a  solution.  The  measured  and  CFIE-computed  radar  cross  sections  of 
a  cone-sphere  with  a  half  angle  of  7°  were  shown  in  [22]  to  disagree  by  several  dB  over  a  large 
range  of  incident  angles,  and  it  remained  uncertain  as  to  whether  the  measured  data,  the 
CFIE  computations,  or  both  contained  large  errors.  Results  for  this  cone-sphere  computed 
with  the  DSEFIE  developed  in  the  present  paper  agreed  very  closely  with  the  CFIE  results 
in  [22],  and  thus  it  was  concluded  that  the  large  discrepancy  between  the  computed  and 
measured  radar  cross  sections  was  due  mainly  to  errors  in  the  measured  data. 

The  importance  of  having  two  independent  surface  integral  equation  formulations  is  also 
underlined  by  the  fact  that  calculating  scattering  with  surface  integral  equations  is  often  far 
from  being  a  simple  “turn  the  crank”  procedure.  As  a  recent  investigation  of  the  convergence 
properties  of  the  CFIE  has  shown  [23],  the  value  of  the  CFIE  combination  parameter  ao 
can  strongly  influence  the  rate  at  which  scattering  calculations  converge  with  increasing 
grid  point  density.  Similar  conclusions  apply  to  our  own  experience  with  the  influence  of 
the  combination  parameter  a  and  the  separation  distance  S  on  results  obtained  with  the 
DSEFIE.  Careful  checking  of  the  results  of  scattering  calculations  obtained  with  surface 
integral  equations  is  thus  very  important,  and  can  be  done  more  easily  and  rigorously  if  two 
independent  formulations  are  available. 

®The  difficulty  with  the  MFIE  and  DSMFIE  near  the  tip  of  a  small-angle  cone-sphere  several  wavelengths 
long  has  been  remedied  by  inserting  the  known  Bessel  function  dependence  of  the  exact  zeroth  order  lon¬ 
gitudinal  tip  current  of  the  infinite  cone  [20,  ch.  5].  There  is  no  guarantee,  however,  that  this  remedy  will 
work  for  any  shape  and  size  of  a  BOR  containing  a  narrow  tip. 
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The  organization  of  the  report  is  as  follows.  Section  2  contains  the  analysis  of  the  DSEFIE 
solution  of  the  BOR  scattering  problem.  It  is  divided  into  several  subsections  beginning  with 
the  statement  of  the  problem  and  definition  of  the  geometry  in  Section  2.1.  The  formulation 
of  the  DSEFIE  in  terms  of  the  vector  and  scalar  potential  functions  A  and  $  is  derived  in 
Section  2.2.  The  solution  of  the  DSEFIE  by  the  Galerkin  form  of  the  method  of  moments  is 
outlined  in  Section  2.3.  Detailed  expressions  for  the  elements  of  the  Z  matrices  that  multiply 
the  column  vectors  of  the  surface  current  expansion  function  coefficients  to  be  determined 
are  derived  in  Section  2.4.  The  Z  matrices  are  treated  by  expressing  them  as  the  sum  of  four 
submatrices.  Two  of  these  submatrices  contain  the  terms  involving  the  vector  potential  A, 
one  each  for  the  observation  point  on  the  original  surface  and  the  dual  surface,  respectively, 
and  two  contain  the  terms  involving  the  scalar  potential  $,  again  one  each  for  the  observation 
point  on  the  original  and  dual  surface.  It  is  noteworthy  that  the  procedure  of  Mautz  and 
Harrington  [11]  we  use  to  transfer  the  differential  operator  on  $  to  the  testing  function  when 
the  observation  point  is  on  the  original  surface  cannot  be  used  when  the  observation  point 
lies  on  the  dual  surface.  In  Section  2.5  we  obtain  detailed  expressions  for  the  elements  of  the 
V  column  vectors  in  the  right-hand  side  of  the  Galerkin  matrix  equation  formulation  of  the 
DSEFIE.  In  Section  2.6  we  obtain  expressions  for  the  currents  induced  on  the  surface  of  a 
BOR  by  a  transverse  electric  (TE)  and  transverse  magnetic  (TM)  linearly  polarized  plane 
wave  in  terms  of  the  solution  to  the  Galerkin  matrix  equation,  and  in  Section  2.7  expressions 
for  the  components  of  the  far  scattered  field  are  derived.  The  analysis  of  Sections  2.4  -  2.7  is 
general  in  the  sense  that  no  explicit  form  is  assumed  for  the  expansion  and  testing  functions 
other  than  the  Fourier  modes  used  to  express  their  azimuthal  dependence.  In  Section  2.8 
we  introduce  the  four-impulse  approximation  to  a  triangle  function  used  to  express  the 
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dependence  of  the  expansion  and  testing  functions  on  the  BOR  generating  curve  parameter 
t.  This  allows  the  integrations  with  respect  to  t  of  the  expansion  and  testing  functions  to  be 
carried  out  in  closed  form.  The  expressions  for  the  elements  of  the  Z  matrices  and  V  vectors 
obtained  in  Sections  2.4  and  2.5,  and  the  expressions  for  the  far  scattered  field  obtained 
in  Section  2.7  are  then  given  as  summations  that  can  be  readily  evaluated  by  computer. 
The  efficiency  of  evaluating  these  summations  can  be  greatly  increased  by  using  a  change  in 
indexing  described  in  Section  2.9.  Section  2.10  discusses  the  choice  of  the  number  of  Fourier 
modes  that  need  to  be  used  in  the  calculations. 

Section  3  contains  the  results  of  calculations  performed  with  a  DSEFIE  computer  pro¬ 
gram  written  to  implement  and  validate  the  analysis  of  Section  2.  Section  3.1  shows  that  the 
DSEFIE  removes  a  spurious  resonance  that  appears  when  the  backscatter  radar  cross  section 
(RCS)  of  a  PEC  sphere  is  calculated  with  the  EFIE.  Section  3.2  demonstrates  the  DSEFIE 
removal  of  a  spurious  resonance  in  the  bistatic  RCS  pattern  of  a  PEC  prolate  spheroid  ob¬ 
tained  with  the  EFIE.  Section  3.3  shows  the  effectiveness  of  the  DSEFIE  in  eliminating  a 
spurious  resonance  in  the  calculation  of  the  bistatic  and  monostatic  patterns  of  a  finite  cylin¬ 
der.  Section  3.4  discusses  calculations  of  the  bistatic  RCS  pattern  of  an  axially  illuminated 
PEC  cone-sphere,  and  the  special  treatment  of  the  dual  surface  needed  when  the  DSEFIE 
is  applied  to  a  cone-sphere  with  a  narrow  tip  angle.  Section  3.5  continues  the  discussion  of 
scattering  by  a  cone-sphere  with  a  narrow  tip  angle  by  treating  monostatic  scattering.  The 
monostatic  case  is  significantly  different  because  of  the  presence  of  the  zeroth  order  Fourier 
mode,  absent  in  axial  illumination.  Section  3.6  discusses  the  use  of  half-triangle  basis  func¬ 
tions  at  the  beginning  and  end  of  the  generating  curve  for  a  BOR.  A  report  summary  is 
given  in  Section  4. 
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2  ANALYSIS 


2.1  Statement  of  Problem 

We  seek  to  determine  the  surface  current  and  the  far  scattered  field  of  a  perfectly  electrically 
conducting  (PEC)  closed  body  of  revolution  (BOR)  excited  by  an  incident  plane  wave.  The 
geometry  of  the  BOR  is  shown  in  Fig.  2.  Circular  cylinder  coordinates  {p,  </>,  z)  are  employed 
with  (p,  z)  denoting  the  corresponding  unit  vectors,  and  with  the  2:  axis  chosen  as  the  axis 
of  revolution.  The  origin  of  the  circular  cylindrical  coordinate  system  lies  on  the  z  axis  but 
does  not  necessarily  coincide  with  the  lower  pole  of  the  BOR  as  in  Fig.  2.  The  coordinates 
(t,  0),with  t  the  pathlength  along  the  generating  curve  of  the  BOR  from  the  lower  pole, 
form  an  orthogonal  curvilinear  system  on  the  surface  S  of  the  BOR;  the  corresponding  unit 
vectors  are  (t,  ^).  Figure  3  shows  the  propagation  vector  =  k  of  the  incident  plane 
wave.  The  propagation  vector  is  assumed  to  lie  in  the  xz  plane  (^  =  0),  with  — k*"^  making 
an  angle  of  0*”^^  with  the  positive  z  axis  and  with  <  0  so  that 

k*"'^  =  -A:  (sin  +  cos  0*”"z) .  (10) 

The  free  space  propagation  constant  is  given  by  k  =  ufc  where  c  is  the  speed  of  light 
and  a;  >  0.  Harmonic  time  dependence  exp(juit)  is  assumed.  Also  shown  in  Fig.  3  are 
the  spherical  polar  angles  of  the  far  field  observation  point  =  (r,  and  the 

/V  fav  ^  fCLT 

associated  unit  vectors  6  and  0  .  For  TM  illumination  the  incident  electric  field  is 

given  by 

E*""  =  fcZoexp(-ik*”‘=  •  r)e*""  (11) 
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Figure  2.  Body  of  revolution  and  coordinate  system. 
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Figure  3.  Plane  wave  scattering  by  a  body  of  revolution. 
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while  for  TE  illumination 


=  fcZoexp(-ik‘""  •  (12) 

In  (11)  and  (12)  r  is  the  vector  from  the  origin  to  any  point  in  space  and  Zq  is  the  intrinsic 
impedance  of  free  space.  The  factor  of  kZo  is  inserted  to  simplify  the  analysis. 


2.2  Derivation  of  the  Dual-Surface  Electric-Field  Integral  Equa¬ 
tion 


To  derive  the  dual-surface  electric-field  integral  equation  (DSEFIE)  we  first  derive  the  or¬ 
dinary  electric- field  integral  equation  (EFIE).  On  the  surface  S  of  the  PEC  BOR  the  total 
tangential  electric  field  vanishes.  The  total  field  is  expressed  as  the  sum  of  the  incident  field 
and  the  scattered  field,  so  that 

n(r)  X  E‘°‘(r)  =  n(r)  x  [E*"^(r)  -H  E*^(r)]  =  0,  r  on  5  (13) 

where  E‘°‘  and  E**^  are  the  total  and  scattered  electric  fields,  respectively,  and  n  is  the  unit 
normal  vector  to  the  surface  S  at  r,  assumed  directed  outward  from  S.  The  scattered  field 
can  be  expressed  in  terms  of  a  vector  potential  A(r,  J)  and  a  scalar  potential  $(r,  J)  by 

E^‘^(r)  =  -ja;A(r,  J)  -  V«>(r,  J)  (14) 


where 


and 


(15) 


(16) 
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Here  <j(r,r')  is  the  free-space  Green’s  function 


_/x  exp(-ifc  I  r  -  r'  I) 

r  j  = - — j - T] - 1 

47r  r  —  r' 


(17) 


r  and  r'  are  the  vectors  to  the  field  and  source  points  respectively,  J(r')  is  the  electric 
current  on  S  to  be  determined,  /xq  and  €o  are  the  permeability  and  permittivity  of  free  space 
respectively,  and  a  is  the  surface  charge  density  given  by 


<'(>■')  =  (18) 

where  the  operator  is  the  surface  divergence  [16,  Appendix  2,  18.].  Combining  (13)  and 
(14)  and  dividing  by  the  free-space  impedance  Zq  =  (/xq/co)^^^  we  obtain 

^  X  [ja;A(r,  J)  -f  V$(r,  J)]  =  x  F-(r),  r  on  5.  (19) 

Equation  (19)  with  A  and  $  given  by  (15)  and  (16)  is  the  potential  form  of  the  EFIE  for 
the  current  on  the  surface  S  of  the  BOR. 

To  obtain  the  DSEFIE  we  note  that  for  points  inside  of  S  the  total  electric  field  vanishes 
so  that 

E"‘=(r)  =  -E'"‘’(r),  r  inside  S.  (20) 


With  (14)-(16),  (20)  can  be  written  as 

iu;A(r,  J)  +  V$(r,  J)  =  E*"‘’(r),  r  inside  S.  (21) 


In  (21)  we  can  now  let  r  lie  on  a  surface  Ss  parallel  to,  and  a  small  distance  <5  >  0  inside, 
the  actual  surface  of  the  BOR.  If  r'^(r)  denotes  the  point  on  the  dual  surface  corresponding 
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to  a  point  r  on  the  original  surface^  then 


ia;A(r‘'(r),  J)  +  Vr<f$(r‘^(r),  J)  =  E'"‘^(r‘'(r)),  r  on  5.  (22) 

Dividing  (22)  by  Zq  and  adding  an(r)  x  (22)  to  (19)  yields  the  DSEFIE 

X  [ju;A(r,  J)  +  jQ!a;A(r‘^(r),  J)  +  Vr$(r,  J)  +  aVH^(r‘'(r),  J)] 

=  ^  X  [E‘”^(r)  +  aE'”^(r‘'(r)l ,  r  on  5.  (23) 

If  the  combination  constant  a  has  an  imaginary  part  and  the  separation  distance  S  is  less 
than  half  a  wavelength  the  DSEFIE  (23)  provides  a  unique  solution  for  J(r)  at  all  real 
frequencies  [13],  unlike  the  conventional  EFIE  (19).  We  have  found  that  choosing  a  =  .25j 
has  worked  well  for  the  BOR’s  we  have  considered. 


2.3  Solution  of  the  DSEFIE  by  the  Galerkin  Form  of  the  Method 
of  Moments 


To  solve  the  DSEFIE  (23)  for  the  surface  current  J  we  employ  the  method  of  moments  and 


let 


^  j 


(24) 


where  and  are  the  expansion  (basis)  functions  of  the  BOR  surface  coordinates  t  and 


0, 


Jni  =  n  =  0,  ±1,  ±2,  •  •  • ,  (25a) 

"^For  smooth  generating  curves  can  be  obtained  from  r  by  setting  r'^(r)  =  r  —  (Jn(r).  For  generating 
curves  with  corners,  however,  a  correspondence  such  as  shown  in  Fig.  4  for  the  generating  curve  of  the  finite 
cylinder  may  be  preferable. 
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Figure  4.  Correspondence  between  points  on  the  body  of  revolution  generating  curve  for  a 
finite  cylinder  and  points  on  the  dual-surface  generating  curve. 


<•  =  n  =  0,  ±1,  ±2,  ■  • . ,  {25b) 

and  are  coefficients  to  be  determined,  and  the  summation  over  the  Fourier  modes 
jg  jy  rpj^g  choice  of  the  functions  fj{t)  is  treated  below  in  Section  2.8,  and 

the  value  of  N  discussed  in  Section  2.10. 

Next,  using  the  Galerkin  method,  the  dot  product  of  (23)  with  each  one  of  a  collection 
of  testing  functions  Wj„,,  defined  by 


WL  =  ifi{t)e-^^^,  m  =  0,  ±1,  ±2, . . . ,  ±N, 

(26a) 

Wii  =  0/i(t)e-^-^  m  =  0,  ±1,  ±2,  •  •  • ,  ±iV 

(266) 

is  integrated  over  S  and  use  is  made  of  the  fact  that  n,  0,  t  form  an 

orthogonal  triad  of 

vectors  to  yield  the  matrix  equations 

+  [z'JM  =  Vi. 

n 

(27a) 

Y.azS.K  +  [ztM)  =  Vi 

(276) 

n 


where  the  Z's  are  square  matrices  whose  elements  are  given  by 

(Zalij  =  /W^,  •  [ji.;A(r,  j;,)  +  jaa.A(r''(r),  j;,) 

+V,«(r,  j;,)  +  oV,.'l.(r‘'(r),  (28) 

Vi.  and  V*  are  column  vectors  whose  i**  elements  are 

C  =  ^  /  WJ. ,  ■  [e‘”  +  QE™'(r''(r))]  dS.  (29) 

and  1^,1^  are  the  column  vectors  of  the  coefficients  of  the  expansion  functions  in  (24).  In 
(28)  and  (29)  p  and  q,  which  can  be  either  t  or  (f),  indicate  the  unit  vectors  of  the  testing 
and  expansion  functions,  respectively. 
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2.4 

Expressions  for  the  Elements  of  the  Z  Matrices 

Detailed  expressions  are  now  obtained  for  the  elements  of  the  Z  matrices  given  by  (28).  We 

write 

+ o‘\zz.t  +  [zai* + 

(30) 

with 

raj  =  ;^  /  WL  ■  jwA(r,  jyrfs, 

s 

(31a) 

\Z2Xt‘  =  ^  /  WL  •  Ju,A(r''(r), 

(316) 

My  =  ^  /  WL  •  V,4.(r, 

Z^J 

(31c) 

(My-  =  Yo 

s 

(31d) 

and  treat  each  of  these  four  terms  in  turn.  The  superscript  d  is  used  here  and  below  to 
indicate  dependence  on  the  dual  surface.  The  superscripts  p  and  q  indicate  the  unit  vectors 
of  the  testing  functions  and  the  expansion  functions,  respectively. 


Equation  (31a)  together  with  (15),  (17),  (25),  and  (26)  becomes 

TT  TT 

— TT  — TT 

In  (32)  and  the  following,  quantities  with  primes  are  associated  with  the  source  point  (ex¬ 
pansion  function)  integration  and  unprimed  quantities  are  associated  with  the  observation 
point  (testing  function)  integration.  The  integrations  with  respect  to  t  and  t'  are  from  0 
(the  lower  pole  of  the  BOR,  see  Fig.  2)  to  the  length  of  the  generating  curve  for  the  BOR, 
the  value  of  t  and  t'  at  the  upper  pole  of  the  BOR.  It  is  easy  to  establish  the  unit  vector 
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relations^ 


and 


t  •  =  sin  v(<)  sinu'(<')  cos{(f>'  —  4>)  +  cosv{t)  cosv'{t'), 

t  •  0  =  —  sin  v{t)  sin((^'  —  (p), 

^  =  sin  v'{t')  sin(<^'  —  (p), 

<p  -  ^  —  cos{(p'  —  <p), 


where  v  and  v'  are  the  angles  measured  positive  clockwise  from  the  positive  z  axis 


t  respectively  (see  Fig.  5).  It  is  also  simple  to  show  that 


r  -  r  h= 


-  2pp'  cos{(p  -  (p')  +  {z-  2')^] 

=  [(P  -  pT  +  (2  -  z'f  +  4pp'sin2(^-J-^)]^^^ 

In  view  of  (33)  and  (34)  the  cp  and  (p'  integrations  in  (32)  are  of  the  form 

TT  TT  TT  TT 

j  d(pe-^'^^  I  d(p'e^^^' f{<p'  -  0)  =  /  d(pe^^^-^^^  j  d(P' f{(p'  -  <p) 

—  IT  — TT  —TT  —TT 

TT 

=  j  d4,’e’'‘*  S(4,’) 

—  TT 

where  the  Kronecker  delta  6nm  equals  0  for  m  /  n  and  equals  1  for  m  =  n,  and 
used  the  27r-periodicity  of  the  ^'-integrand  to  set  0  =  0.  Hence 

TT 

\ZZa\t,  =  \Z”]tj=i2irk  jdtpMt)  J  dt'p'/,(t')  J  d4,'{t,  ci’)e”'*'G(R) 


^Equations  (33a-d)  follow  directly  from  noting  that 


t  =  cos  v{t)z  +  sin  v{t)p  =  sin  v{t)  cos  0x  +  sin  u(^)  sin  </>y  +  cos  z;(f)z, 


^  =  —sin4>x  +  cos  0y, 


and  similarly  for  t  and  (f) . 


(зза) 

(ззб) 
(33c) 

(33d) 
to  t  and 

(34) 


(35) 
we  have 

(36) 
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Figure  5.  Defining  geometry  for  the  angle  v  in  the  pz  plane. 
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where,  from  (33)  and  (34)  with  (/>  =  0, 


t  •  =  sin  vit)  sin  v'it')  cos  (j)'  +  cos  vit)  cos  v'it'), 

(37a) 

t  •  0  =  —  sin  vit)  sin  (f) , 

(37(.) 

(})  -t  =  sin  v'it')  sin  4> , 

(37c) 

4>  •  (f>  =  cos  (f)  , 

(37<i) 

r  1  1/2 

i?  =1  r  -  r'  1=  (p  -  p')2  +  (z  -  +  4pp'sin2(|)  . 

(38) 

Then 

=;27rA:  J  dtpfi{t)  j dt'p'fj{t') 

TT 

•  j  d4>G{R)e^'^'^' cos (})'  ■\-cosv{t)  cos  v'{t')],  (39o) 

—  TT 

7r 

=  j2'Kk  j  dtpfi{t)  j  dt'p' fj{t')  sin  v'{t')  J  d<j>'G{R)e^’^'^' sin  <f)',  (396) 

—  TT 
TT 

[Z^]^  =  —j2‘Kk  j  dtp fi{t)  sin  v{t)  j  dt'p'fj{t')  J  d0'G'(/?)e^"'^' sin  </)',  (39c) 


=  >27rA:  I  dtpfiit)  I  dt'p'fjit')  J  #'G(i?)e^"^'  cos  (f>'.  (39d) 

—  TT 

Making  use  of  the  symmetries  in  the  <f>'  integrations  in  (39),  we  obtain 

[Z“]J  =  jk'^  I  dtpfiit)  J  dt' p'f j{t')[sinvit)  sin  v'{t')G2, nip,  p\  z  -  z') 

+  cost;(t)  cosr;'(t')Gi,n(p, -2  -  20]» 

(40a) 
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and 


where 


and 


with 


=  -k^  J  dtpfiit)  j  dt' p' sin  v'{t')G3,n{p,  p' ,  z  -  z'). 

(406) 

l^n  ]ij  =k^  1  dtp fi{t)  sin  v{t)  j  dt'p' fj{t')G3,n{p,p',z  -  z'). 

(40c) 

[Zt^]fj=^jk^l  dtpfiit)  1  dt'p'fj{t')G2,n{P,P',z-z'), 

(40d) 

TT 

Gi,n{p,p',  Z-  z')  =  j  GoiR)  cos{n<f)')d<l)' , 

0 

(41a) 

TT 

<J2,n(p,  p',  z-  z')  =  j  Go{R)  cos{n<f>')  cos  (f)'d(f)', 

0 

(416) 

TT 

G3,„(p,  p',  z-  z')  =  j  Go{R)  sin{n<j)')  sin 

0 

(41c) 

^  /m  exp{-jkR) 

^  kR 

(41<i) 

and  R  given  by  (38).  In  the  arguments  of  the  Gk,m  k  =  1,2,3,  p  and  z  are  functions  of  the 
observation  point  coordinate  t,  and  p'  and  z'  are  functions  of  the  source  point  coordinate 
t'.  Equations  (40a-d),  the  detailed  expressions  for  (31a-d),  are  in  terms  of  the  t  dependence, 
fi{t),  of  the  expansion  and  testing  functions  discussed  below  in  Section  2.8.  Referring  to 
(38)  it  can  be  seen  that  R  =  0  when  p'  =  p,  z'  —  z,  and  =  0,  so  that  Gi,n(p,  p',  -  z')  and 
G2,n{p,  p',  z  —  z')  become  singular.  This  problem  will  also  be  treated  in  Section  2.8. 

Turning  next  to  (31b)  it  can  be  seen  that  the  treatment  of  (31a)  carries  over  if  = 
{p‘^,(j),  z‘^)  is  substituted  for  r  in  (32)  and  (34).  (Note  that  the  value  of  0  is  the  same  for 
corresponding  points  r  and  r‘^(r),  and  that  the  testing  function  integrations  are  performed 
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on  the  surface  of  the  BOR.)  Hence  letting 


=1  -  r'  I,  (42) 

[Z“]§‘'  =  jk'^  j dtpfiit)  j dt'p'fj{t')[sm  v{t)  sin  t;'(t')G2,n(p^ p\  -  z') 

+  COSv{t)cOSv'{t')Gi^n{p'^,p\z'^  -  z')]. 

(43a) 

[Zf  1^5''  =  -k'^  j  dtpfiit)  j  dt'p'fjit')  sin  v'{t')Gz,n{p‘^ ,  p\  z‘^  -  z'),  (436) 

=  k'^  j  dtpfi{t)  sin  v{t)  j  dt' p'  fj{t')Gz,n{p'^,  p',  z^  -  z'),  (43c) 

and 

[Z**]'^"  =  jk-‘ I  dtpMt)  W 

where  the  Gk,n,k  =  1,2,3,  are  defined  as  in  (41a-d)  with  R  replaced  by  R^.  It  should  be 
noted  that  unlike  R  defined  by  (38),  R'^  defined  by  (42)  cannot  equal  zero  and  so  there  are  no 
singularities  of  the  integrands  of  Gi,n(p‘^,  p',  z^  -  z'),  G2,n(p‘^,  p',  z'^  -  z'),  and  G^^nip^,  p\  z'^  - 

Next,  addressing  (31c),  we  first  show  that  [11] 

/  WL  ■  V,4'(r.  =  -| (Vj  ■  WL)4>(r,  j;,)<iS,  (44) 

s  s 

where  Vs-  is  the  surface  divergence  operator  [16,  Appendix  2,  18.].  The  importance  of  (44) 
is  that  the  integrable  singularity  of  |  r  -  r'  [“'  in  the  Green’s  function  G(r,  r')  given  by  (17) 
is  thereby  not  converted  to  a  singularity  of  higher  order  |  r  — r'  by  the  operator  Vr  acting 
on  4>  in  the  left-hand  side  of  (44).  To  prove  (44)  we  start  with  the  relationship 

I  Vs  ■  (4’W^,)(iS  =  0  (45) 

s 
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which  holds  because  S  is  closed  and  W^j  is  tangential  to  S  [16,  Appendix  2,  42.].  Since 


Vs  ■  =  WL  •  Vs4  +  (Vs  ■  WL)$,  (46) 


we  have 


I WL  •  VsidS  =  -l(Vs-  W'„,)idS  (47) 

5  S 

and  (44)  follows  from  noting  that  in  (47)  Vs^  can  be  replaced  by  V<I»  because  W^j  is 
tangential  to  S.  Equation  (31c)  together  with  (44), (16), (17),  and  (18)  then  becomes 


Now 


and 


Udsf  dffiVs  ■  WL)(V's  •  j;,)G(r,  r'). 

(48) 

iv  J  J 

s  s 

V's-Ji,  =  — 

(49a) 

(496) 

Vs-Wi.,  =  i|H(()l<J’‘*', 

(49c) 

Vs-Wt  =  -^/i(i)«4”*'. 

P 

(49d) 

Hence 


[ZLlti  =  -{!  I  I  dt'--^[p'm]  I  d<l>'p'e^^‘f>’G{\  r  -  r'  |), 


(50a) 


[ZL]ti  =  dt^-^  I  d<j>pe-^'^^j  dt'--[p'm]  I  r  -  r'  |),  (506) 


—  TV 
TV 


=  11  I  d4,pe-’^*  j  d(fi^  J  dd.'p'ei’l'ai]  r  -  r'  |),  (50c) 
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and 


iztfjt,  =  /  d4,pe-’"‘*fdt'li^  j  #VV"*'G(|  r  -  r'  |).  (50rf) 

—  TT  — TT 

Just  as  for  the  </»  and  integrations  in  (50)  are  treated  using  (35)  and  the  symmetries 

of  the  resulting  <f)'  integrations  then  yield 

[^n]ij  =  -j  P'^  2  -  ^')>  (51a) 

[Zt%  =  -nl  dt^  j  dt'-^[p'm]Gr,n{p.  p',  z  -  z'),  (516) 

[Zn%  =  n  j  dt^[pfi(t)]  j  nip,  p',  z  -  z'),  (51c) 

and 

|zri?  =  J  dt^  j df^^GUP,  p',  2  (sw) 

Finally  we  turn  our  attention  to  (31d).  Here  it  is  not  possible  to  use  the  device  of  (44) 
to  transfer  the  differential  operator  on  <F  to  the  testing  function  W^j  because  the  gradient 
operates  with  respect  to  points  on  the  dual  surface.  Hence,  using  (16), (17),  and  (18),  we 
obtain 


(52) 


SO  that 


dS'lV’s  ■  ■  (f  -  r').  (53) 


1  +  JA:  I  -  r'  I 

uo  J  UJ  [\  s  ■  Onj\r 

s  s 

It  is  simple  to  establish  the  relations  (see  the  footnote  regarding  (33a-d)) 


t  •  (r*^  -  r')  =  p^  sin  v{t)  —  p'  sin  v{t)  cos{<f>  -  0)  +  {z’^  —  z')  cos  vit)  (54o) 
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Hence  (53)  together  with  (54)  and  (49a, b)  yields 


—  TT  — TT 

p**  sin  v{t)  —  p'  sin  v{t)  cos{<f>'  —  <^)  +  {z^  —  z')  cos  v{tj^  H{R^), 


(55a) 


[Ztn%  =  I  j  dtfiit)  I  d<f>pe-^^  j  dt'--[p'fj{t')]  I  d<i>'p'^e^^^'sm{<f>'-<l>)H{R%  (555) 


=  ll  dtfiit)  I  d<l>pe-i^‘>>  I  df-m  I  #'p'e^”"' 

— ’TT  — TT 

[p*^  sin  D(t)  —  p'  sin  v{t)  cos(^'  —  ^)  +  {z‘‘'  —  z')  cos  u(t)j  H{R^), 


(55c) 

and 

TT  -  TT 

[Zttn]^  =  -^/  di/i(0  /  d<j>pe-i^^  j  I  d<j>’p'h^^’»'sm{<f>'  -  cl>)H{R%  (55d) 

—  TT  “TT 

where  we  have  let 

Tj(p<i\ _  ^  jkR  jkuA 

with 

=1  r**  -  r'  I  .  (57) 

Once  again  the  (j)  and  (f)'  integrations  are  treated  using  (35),  and  the  symmetries  of  the 
resulting  (^'  integrations  then  yield 


[Zn]ij  =  j  I  dt' -^,[p' fj{t')] 

•  {[p‘'sin?;(^)  +  (2*^  -  2')cost;(<)]  ,  p\  z'^  -  z')  -  p'  smv{t)H2^n{p‘^,p\z^  -  2')}  , 


(58a) 


[Ztt  =  -k'"  I dtpfiit)  j 


[Zn\t  =  dtpfi{t)  j  i 


jp'fjit') 


{[p‘'sinv(f)  +  {z"^  -  z')  cosu(t)]  //i,„(p'',p', 2*^  -  z')  -  p' smv{t)H2,n{p‘^,p\z’^  -  2')}  , 


In  (58) 


where 


(58c) 

]tj  =  -jnk'^  J  dtpfiit)  J  dt'p'fj{t')H3,n{p‘^,p',z^  -  z'). 

(58d) 

TT 

R,n{p‘',p',z^-z')  =  j  Ho{R)cos{n<l>’)d(f>\ 

0 

(59a) 

TT 

H2,n{p‘^,  p\  z.'^  -  z')  =  J  Ho{R)  cos{n<p')  cos  (j)'d<t), 

0 

(596) 

TT 

Hz,n{p^,  p',  -  z')  =  j  Ho{R)  sm{n(j)’)  sin  (j)'d<t>', 

0 

(59c) 

{kRdf 

(59d) 

and  is  given  by  (57). 

In  view  of  the  fact  that  we  have  demonstrated  that 


[^mnh  =  0,  m  /  n 
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the  set  of  matrix  equations  (27)  for  determining  the  coefficient  vectors  IJj  and  reduce  to 

,  n  =  0,  ±1,  ±2,  •  •  • ,  ±A^.  (61) 

It  is  also  easy  to  see  from  (30)  and  (40),  (43),  (51),  and  (58)  that 


- 1 

s 

. . i 

I‘ 

II 

I - 

C 

> 

1 _ 

_ i 

I 

1 

HH 

i _ 

1 

> 

II 

(62a) 

1 

II 

(626) 

1 

II 

(62c) 

2.5  Expressions  for  the  Elements  of  the  V  Vectors 

Now  that  we  have  obtained  detailed  expressions  for  the  elements  of  the  Z  matrices  we  turn 
to  the  elements  of  the  column  vectors  and  in  (27)  given  by  (29).  Let 

VS'=;^/w;,-E‘”«(r)dS  (63) 

where  p  is  either  t  or  (f),  indicating  the  unit  vectors  of  the  testing  functions  given  by  (26), 
and  q  is  either  6  or  <l>,  indicating  the  polarization  of  the  incident  electric  field  given  by  (11) 
or  (12).  The  Fourier  mode  index  m  has  been  replaced  by  n  since  the  original  distinction 
between  the  Fourier  mode  index  n  of  the  expansion  functions  for  the  current  and  the  Fourier 
mode  index  m  of  the  testing  functions  is  no  longer  needed  in  view  of  the  orthogonality 
relation  (35)  which  has  been  used  in  deriving  expressions  in  Section  2.4  for  all  the  elements 
of  the  Z  matrices.  Since  the  expressions  that  will  be  obtained  for  V^i  given  by  (63)  will 
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carry  over  immediately  for  needed  in  (29),  there  is  no  need  to  treat  these  two  cases 

separately.  It  is  simple  to  establish  the  unit  vector  relations  (see  the  footnote  to  (33a-d)) 


i  d  =  cos  0'"^^  sin  v{t)  cos(</>  -  -  sin  0*"^^  cos  u(t), 


^■0  =-cos9^^'^sm{<f>-f^^), 


^  ^  ITIC 

t  •  ^  =  sin  v{t)  sin((^  - 


and 


<i>  •  ^  =  COs{(j)  — 


Also,  using  (10), 


jjtnc  .  j.  _  _/j.^sjn  +  cos  ^‘"‘^z)  •  (p  cos  0x  +  p  sin  (f>y  +  zz) 


=  —kp  sin  cos  (j)  —  kz  cos  9^' 


Hence,  recalling  that  it  has  been  assumed  that  =  0, 


27r 

V;‘f  =  k  j  dtp/i(t)e^^""°*®‘'‘'[cos0*”"sinu(t)  j 


and 


-sin^^”"cosu(t)  I  rf0e'(^psin(?*-cos0-n^) 

0 

27r 

=  -cos6>‘""A:  I J  d0sin 

0 

27r 

V^f  =  k  j  d<p/i(t)e^''^“®^‘'"sinu(t)  J  #sin<^eJ(*^^®‘"®’ 


0'^^  cos  (p—n(p) 


'  cos  (t>-n(p) 


27r 

=  k  J  dtp J  dcf)  cos 


sin  cos  4>—n<l>) 


(64a) 

(646) 

(64c) 

(64d) 


(65) 


(66o) 

(666) 

(66c) 

(66d) 
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Using  the  integral  representation  for  the  Bessel  function 


J^{x)  =  4-  / 

ZTT  J 
0 

it  is  simple  to  obtain  the  integral  representations 

27r 

I  =  -irflJn+iix)  +  J„_i(a:)] 


(67) 


(68a) 


and 


2w 

j  -  7rj"+Vn+l(^)  -  Jn-i{x)] 


(686) 


so  that 


Vni  =  j  dtpfi{t)  |^cos0'"‘^sinn(t)(J„+i  -  Jn-i)  +  2jsin^*"‘^cosi>(t)J„] 


kz  COS  6'^^ 


and 


where  we  have  let 


(69a) 

Kf  =  TTfk  1  dtpfiit)  cos  e^^Vn+i  + 

(696) 

K?  =  -T^fk  j  dtpfiit)  s\nv{t)(Jn-\.i  +  J„_i)e’*^‘'“'’‘"', 

(69c) 

Ki?  =  1  dtp/i(t)(J„+i  - 

(69d) 

Jn  =  Jnikpsme^^^). 

(69e) 

Using  the  Bessel  function  relation 


j_„  =  i-iruz) 


(70) 


it  is  easy  to  show  from  (69)  that 


yte  _  yte 
^  —m  ^  m  5 


(71a) 
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(716) 


^ —ni  ^ni  ’ 

V-i,  =  (71c) 

and 

Kf?.  =  V„**.  (71rf) 

The  expressions  corresponding  to  (69)  for  V^i  when  in  (63)  is  evaluated  at  r‘^(r)  on 
the  dual  surface,  needed  for  calculating  (29),  are  obtained  from  (69)  simply  by  substituting 
for  p  in  (69e)  only  and  for  z  in  (69a-d). 


2.6  Calculation  of  the  Current 


To  calculate  the  current  on  the  BOR  surface  we  refer  to  (24)  and  (25)  and  write 

J’((,0)=  l;  e<"*i;[i'1/i(i)t  +  /S/i(<)0]  (72) 

n—  —  N  i 

where  q  =  6  (TM)  or  <f)  (TE)  indicates  the  polarization  of  the  incident  electric  field  given 
by  (11)  or  (12),  and  are  the  elements  of  the  coefficient  vectors  obtained  as 


the  solution  of  the  matrix  equation  (61) 


Kl 

_ 1 

1 

HH 

_ 1 

\rt<i 

’  n 

_ 1 

1 

109 

,n  =  0,±l,±2,---,±;V. 


(73) 


Letting  /  =  [fi{t)f2{t)  •  •  •  fN{t)]  be  the  row  vector  of  the  expansion  functions  fi{t)  (we  use 
the  tilde  to  indicate  the  transpose  of  a  column  vector),  (72)  can  be  written  as 

J«((,«=  ■£  [/I?t  + /IJ'0] .  (74) 

n=—N 

From  (73),  the  relations  (62)  for  the  matrices  in  terms  of  the  matrices,  and  the 
relations  (71)  for  the  VLn  vectors  in  terms  of  the  Vn  vectors,  it  is  simple  to  derive  the 
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relations 


TtO  _  jtO 

(75o) 

J<l>^  ^  T(l>0 

n  “  5 

(756) 

HH 

II 

1 

1— 1 

(75c) 

and 

T</>^  _  T(t>0 

±-n  ”  ‘*'71  * 

(75<i) 

(Note  that  (75b)  and  (75c)  imply  that  Iq^  =  0  and  =  0.)  Substituting  (75)  in  (74)  then 

yields  the  expressions  for  the  currents  induced  on  the  surface  of  the  BOR  by  a  TM  and  TE 

linearly  polarized  incident  plane  wave,  respectively: 

J®(t,  (p)  =  fll^i  + 2  [(/In )  cos  n4>i  +  j  (/if)  sin  ncp^ 

71  =  1 

J 

(76a) 

and 

J*(i,  0)  =  +  2  f;  [(/X*)  sin  nH  +  iht*)  cos 

71=1 

• 

(765) 

2.7  Calculation  of  the  Far  Scattered  Field 

In  the  far-field  region  the  0  and  (p  components  of  the  scattered  electric  field  can 

from  the  vector  potential  A  given  by  (15)  [24,  p.  281] 

be  obtained 

E-(r)  ^^°°-ja;A(r,J). 

(77) 

Then  with  (15)  and  (17), 

Aj  ji*"  1 

£.c«(r)  'i” -jkZo  j  [J’(r')  ■ 

(78) 

5 
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where  p  :  6  or  <f>,  indicates  the  far  scattered  field  component;  p  :  or  is  the 

corresponding  unit  vector  (see  Fig.  3),  q  :  6  or  (j>,  indicates  the  polarization,  9  (TM)  or 
(p  (TE)  of  the  incident  electric  field;  J’  is  given  by  (74),  r  =  (r,  0^“’’)  is  the  far- field 

observation  point,  and  r'  =  {p',(j)',z')  is  the  source  point  on  the  BOR  surface.  In  the  far 
field,  I  r  —  r'  I  in  the  phase  of  the  Green’s  function  can  be  approximated  by 

I  r  —  r'  f  r  —  p'  sin  6^°’’'  cos(<f>'  —  —  z'  cos  9^“’’  (79) 


so  that 


£;*c,P9(r)  ™  -jkZo^^  j [J’(r')  • 


The  relations  (64)  can  be  used  directly  to  obtain 


^far  _  cos{(j)'  —  —  sin  0-^“’'  COS  u(t')  (81a) 

^  ^far  _  ^  {81b) 

i'  •  =  sin  v{t')  sin(0'  —  <^^“’’),  (81c) 

and 

^  =  cos{(p'  —  0-^“’^).  (81d) 

Letting  the  elements  of  the  vectors  and  be  defined  by 


=  kj  fi{t'){t  ■  (82a) 


and 
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with  p  =  6  or  <j),  (80)  together  with  (74)  can  be  written  as 


E‘‘C’P^{r) 


-jkZoe-^'^'^  J:L 


Airkr 


Y,  (rW  +  R?!*')  (83) 


n=-N 


where,  using  (82)  and  the  27r  periodicity  of  the  0'  integration 


2n 

RZ  =  k  I  dtV7i(t')e''="'“"‘’^“'[cos0^“'-sinT;(t')  J 


-  sin  cos  v{t')  I  ^  (84a) 

0 

27r 

Rffi  =  -cos0^“’-A:  I  (846) 

0 

27r 

R^i^kJ  dtV/i(t')e^''"'^“‘'^“'sinu(t')  j  (i(/)'sin0V(''^'“"‘’^‘‘'^°®^'+"^'),  (84c) 


and 


27r 


(84d) 


0 

Comparing  (84a-d)  with  (66a-d)  we  observe  that  the  corresponding  expressions  are  the  same 
if  0*"*=  of  (66)  is  replaced  by  9^°'^  and  —n  of  (66)  is  replaced  by  n.  Then  (69)  and  the  relations 
(71)  can  be  used  to  yield 


and 


3  9^°''^  sin  v{t'){Jn+l  —  Jn-l)  +  “^j  siu  cos  v{t')Jn^ 

gjfer'cose/®'' 

(85a) 

h  f  *V7i(t')cos«^“'{J„+,  + 

(856) 

/  dt'p'Mt')  sin  v(0(J.+i  +  4-, 

(85c) 

,’■+7  Jdt'p'f,{t'){J„^,  - 

(85d) 

35 


where 


Jn  =  Jn{kp' (85e) 
The  relations  (71)  apply  to  the  as  well,  so  that 


ptO  _  ptO 

(86a) 

p4>0  _  _  p4>0 

^~ni  ^ni  ? 

(866) 

pt4>  _  _  pt4> 

(86c) 

and 

e 

0:^ 

II 

1 

(86d) 

Then,  substituting  (86)  and  (75) 

in  (83)  it  is  straightforward  to  obtain 

psc,e0(-\  jkZoe 

^  i„kr 

Ri*l“  +  2  f;  (r"I“  +  )  cos{n^") 

n~\ 

,  (87a) 

r^oo  sin(n(^^“’-), 

oo 

E-.«»(r)  'It”  (R'X  +  R*"ir)  sin(n^/»'), 

n=l 

(87c) 

and 

^kZoe 

'  '  47rA:r 

N 

Rt^lf  +  2J2  (RifiJf  +  cos(nfli^“0 

n—\ 

(87d) 

The  radar  cross  section  a  is  defined  as 

( 

1  Rsc,pg  12 

[7*’^  =  lim  47rr^7— - — 

r^oo  ^inc,q  2 

(88) 

where  p  :  9  or  (j),  denotes  the  component  of  the  far  scattered  field,  and  q  :  6  (TM)  or  (j)  (TE), 
indicates  the  polarization  of  the  incident  electric  field. 
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2.8  Choice  of  the  expansion  and  testing  functions 


The  detailed  expressions  we  have  obtained  for  the  elements  of  the  Z  matrices  and  the  V 
and  R  vectors  all  involve  the  functions  fi{t)  that  give  the  t  dependence  of  the  expansion 
and  testing  functions  in  (25)  and  (26).  The  choice  of  the  fi{t)  has  been  left  open  so  far 
and  must  now  be  addressed.  Following  Mautz  and  Harrington  [11]  we  use  a  four-impulse 
approximation  to  a  triangle  function  developed  as  follows. 

The  generating  curve  of  the  BOR  is  parametrized  in  terms  of  t,  the  distance  along  the 
curve  measured  from  the  lower  pole  of  the  BOR.  For  each  value  of  t,  the  corresponding  point 
on  the  generating  curve  is  given  by  \p{t),  2:(t)].  A  set  of  P  points  pl,P2,”'i  Pp)  P  odd,  P  >  5, 
is  chosen  to  discretize  the  generating  curve  with  p\  =  (pj,  2*)  the  lower  pole  of  the  generating 
curve  corresponding  to  t  =  0,  Pp  =  {p*p,  Zp)  the  upper  pole  of  the  generating  curve,  and  with 
ti  >  tu-  The  generating  curve  is  approximated  by  straight  line  segments  between  adjacent 
points.  The  midpoints  of  the  approximating  straight  line  segments  are  given  by 


(ft.  ^i)  =  ^^^2^)  .•=  1, 2,  1 


and  the  length  of  the  straight  line  segment  denoted  by 


di  = 


21 1/2 


(89) 


(90) 


For  calculation  purposes  the  discretized  generating  curve  then  completely  replaces  the  origi¬ 
nal  generating  curve  and  the  parameter  t  now  becomes  the  length  along  the  discretized  curve 
from  the  lower  pole  instead  of  the  length  along  the  original  generating  curve  from  the  lower 
pole.  Thus,  for  example,  <(p*)  =  0,  /(p^)  =  di,  t{pl)  =  di  +  ^2)  etc.  The  trigonometric 
functions  sint;(t)  and  cosv{t)  appearing  in  (33),  (40),  (58),  (69),  etc.,  with  v  the  angle  mea- 
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sured  positive  clockwise  from  the  positive  z  axis  to  the  tangent  to  the  generating  curve  in 
the  direction  of  increasing  i  (see  Fig.  5),  are  then  replaced  by  the  discretized  values 


sin  Vi  = 


Pi+i  -  Pi 


(91a) 


and 


cos  Vi  = 


^i+i  -  < 


di 


(91(.) 


A  triangle  function  Ti{t)  spanning  four  segments  of  the  discretized  generating  curve  can  be 


defined  as  in  Fig.  6  where 


0, 


n{t)  =  < 


i-4-1 


d2i-,+d2i’  '2i-l  i  '  S  l2i+, 


(92) 


^2i+3  ^  <  f  <  f* 


[  d2i^\  +  d2i^2  ' 


2i-i-3 


The  first  triangle  function  is  Ti{t)  and  the  last  triangle  function  is  t^p-3)/2-  The  triangle 
functions  overlap  as  shown  in  Fig.  7. 

We  now  let  the  functions  fi{t)  be  defined  by 


kpfiit)  =ri{t). 


(93) 


Then  ^[p/i(t)],  also  required  for  calculating  the  elements  of  the  Z  matrices  and  the  V  and 


R  vectors,  is  given  by 


itwm  =  Irm 


(94) 


where 


rlit)  = 


0,  t<t2i-l^t>t2i+3 

^2.-1  <  ^  <  ^2... 


(95) 


I  d2i-M  +  d2i+2 


.  ^2i+l<^<<2i 


2i+3 
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Figure  7.  Overlapping  of  the  triangle  functions  Tj(t)  and  r<+i(t). 
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With  the  functions  fi{t)  chosen  as  above  there  still  remains  the  problem  of  how  the  t 
and  t'  integrations  in  the  expressions  for  the  elements  of  the  Z  matrices  and  the  V  and 
R  vectors  are  to  be  performed,  since  the  use  of  direct  numerical  integration  to  evaluate 
these  expressions  is  clearly  not  practical.  To  perform  these  integrations  we  therefore  use  a 
four-impulse  approximation  to  Ti{t)  (see  Fig.  6)  letting 

4 

kpfi{t)  «  X!  Tp+4i-4S{t  -  tp+2i-2)  (96) 

P=1 


where  the  T's  are  equal  to  the  value  of  r(t)  at  the  segment  midpoints  multiplied  by  the 
width  of  the  segments: 


rp  _  ^t-1 

2{d2i-i  +  d2i)  ’ 

(97o) 

rr  +  \d2i)d2i 

=  ,  ,  ,  > 

+  «2i 

(976) 

rr,  (d2t+2  +  \d2i-i-l)d2i+l 

J4i-1  -  j  ,  j  ’ 

«2i+l  +  “2i+2 

(97c) 

rp  _  ^i+2 

2{d2i+l  +  d2i+2)  ' 

(97d) 

and  6{x)  is  the  Dirac  delta  function.  Similarly,  a  four-impulse  approximation  is 

also  employed 

for  the  ^[pfi{t)]  using  the  derivative  of  r(<)  at  the  segment  midpoints  multiplied  by  1/k  and 

the  width  of  the  segments.  Thus 

(98) 

where 

rp/  _  ^2t— 1 

Cl2i~l  “r  «2i 

(99a) 

rp,  _  d2i 

^  d2t-i  +  d2t’ 

(996) 

rp,  _  d2i+l 

^  ^21+1  -t-  d2i+2  ’ 

(99c) 
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and 


rp!  _  ^21-^2 

d2i+i  +  do 


(99d) 


^2t+l  -r  u2t+2 

As  an  example  of  the  implementation  of  the  four-impulse  triangle  function  approximation, 
let  us  take  (40d) 


[Zt%  =  /  dtpfi{t)  j  dt' p' fj{t')G2, nip,  p',  Z  -  z') 


with  G2,n{p,p',z  —  z)  given  by  (41b).  Substitution  of  (96)  in  (100)  then  gives 


4  4 


lZiX=3Y,Y.Tp’TfG2Ap,;Pf,^,'  -  ^f) 


p=l  q=l 


where 

p'  =  p  -f-  4i  —  4, 
q'  =  q  +  4j  -  4, 
i'  =  p  +  2i  —  2, 
and 

f  =  q  +  2j  -  2. 


In  the  integral  (41b)  that  defines  G2,n,  R  is  given  by 


R  = 


{pi>  -  pj'f  +  (2i-  -  Zj>f  +  4pj'pj/  sin^ 


(100) 


(101) 

(102a) 

(1026) 

(102c) 

(102d) 

(103) 


In  (103),  (pt',2i')  and  iPj',Zji)  are  the  respective  midpoints  of  the  and  segments  of 
the  discretized  generating  curve. 

Once  the  integrations  over  t  and  t'  have  been  replaced  by  summations  as  in  (101),  we 
introduce  the  following  changes  in  the  summation  indexing  in  preparation  for  the  increase 
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in  efficiency  of  the  resulting  calculations  discussed  in  Section  2.9.  First,  let  the  summations 
be  represented  schematically  in  general  as 

=  (104) 

p=lq=l 

where  Qij  is  the  quantity  to  be  evaluated,  is  the  summand  apart  from  the  triangle 

functions,  and  p',  q',  i',j'  are  given  by  (102).  (Although  we  have  used  T  given  by  (97)  in  (104), 
either  or  both  of  Tpf  and  Tq'  can  be  replaced  by  a  derivative  T!^  or  given  by  (99)  without 
in  any  way  affecting  the  argument.)  Next,  we  change  the  indices  so  that  the  summations  are 
over  the  indices  of  F  by  letting 


r  =  i'=p  +  2i  — 2 

(105o) 

so  that 

s  =  j'  =  q  +  2j -2 

(1056) 

p  =  r  —  2i  + 2 

(106o) 

q  =  s  —  2j  +  2 

(1066) 

and 

p'  =  r  +  2i  —  2 

(107a) 

q*  =  s  +  2j  —  2. 

(1076) 

Then  successive 

values  of  p  from  1  to  4  correspond  to  successive 

values  of  r  from  2i  —  1  to 

2i  +  2,  and  successive  values  of  q  from  1  to  4  correspond  to  successive  values  of  s  from  2j  —  1 

to  2j  +  2,  and 

2i+2  2j+2 

<3«=  E  E  Tr.T,F(r,s) 

r=2i-l  s=2j-l 

(108) 
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where 


r'  =  r  +  2i  —  2 


(109a) 


s'  =  s  +  2j  —  2. 


(1096) 


Finally,  we  change  i,j  to  p,  q\  r,  s  to  z,  j;  and  r',  s'  to  i',j'  so  that  (108)  becomes 


2p+2  29+2 

0p.=  E  E  T,T,.F(i,j] 

i=2p-l  3=2q-\ 


(110) 


where 


i'  =  i  +  2p  —  2 


(111a) 


/  =  i  +  2g  -  2. 


(1116) 


Substituting  the  expressions  (96)  and  (98)  for  the  expansion  functions  and  their  deriva¬ 
tives  in  (40), (43), (51),  and  (58),  and  performing  the  change  in  summation  indexing  of 
(104)-(111),  the  elements  of  the  Z  matrices  now  become  the  following: 

2p+2  29+2 

[^n  1m  =  jjl  £  [si"  si"  VjG2,n{pi,  Pj,  Zi  -  Zj)  +  COS  Vi  COS  VjGi,n{Pi,  Pj,  Zi  -  Zj)]  , 

i=2p-l  j=2^-l 

(112a) 


2p+2  29+2 

Ti'Tji  sin  VjGj^riiPi^  Pji  Zi  —  Zj), 

i=2p-l  i=2q-\ 

2p-)-2  29+2 

\Z'^]fj  =  H  iz  Tj>  sin  ViG3,n {pi ,  Pi ,Zi-  Zj), 
i=2p—l  j=2q—l 


(1126) 


(112c) 


2p+2  29+2 

[ZtX,  =j  E  E  T,Tj,G2,n{p^,  Pi,  -  ^i); 

i=2p-l  j=2q-\ 


(112d) 
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2p+2  29+2 

[^n  =  i  S  £  ,  P;,  2:f  -  Zj)  +  COS  Vi  COS  t;jG'i,„(pf,  Pj,  zf  -  Zj)\  , 

i=2p— 1  j=2q—l 

(113o) 

=-  E  E  Ti'Tj>  sin V jGz, n{pf,  Pj,  zf  -  Zj), 

i=:2p-l  j=2q- 
2p+2  2g+2 


2p+2  2(?+2 

]A^  _  rp  rr  ..  n  ( ^ 

Jp9 

z=:2p-l  j=2g-l 


K1$‘=  E  E  T,Ti.smViCsApf.Pi'4-^l), 

i=2p-l  j=2g-l 


and 


A-i 

n  Jp9 


2p+2  2g+2 


JP9  j^2 


n 


and 


«!$<*  _ 


}pq  ^2  ^  /I . 

^  i=2p-^l  jzz2q-l  n 


2  2p+2  2^+2  rp  ji 

-  ~3-jZ  YL  E  )  > 

''  i=2p-li=29-l  PjPJ 


2p+2  2^4-2 

-i  E  E 

i=2p— 1  j=2qf— 1 


(1136) 

(113c) 


E  E  '^i''^j'^2,n{Pi,Pj,zf  -  Zj)] 

i=2p-l  j=2q~l 

(113d) 

2pH-2  2^-|~2 

r  T^tTjtGi^niPii  Pjy  ^  ^j)y 

i=2p—l  j=2q—l 

(114a) 

2p+2  2qi+2 

j  E  E  G'i,n(pnPj,^t  2j), 

'  t=2p-l  j=29-l 

(1146) 

2p-f  2  2^+2  rpf  rp  ^ 

E  E  ‘  /  ^l,n(Pi,Pj,2i  2:j), 

(114c) 

(114d) 


•  { [pf  sin  Vi  +  (zf  -  Zj)  cos  Vi]  /^i,„(pf,  Pj,  zf  -  Zj)  -  p,  sin  ViH2,n{pt  Pj,  4  “  ^j)} 

(115a) 
(1156) 


2p+2  29+2 

i-2p-l  j=2q-l 


[Z^lp, 


2^2  ^2  j,,j, 

"EE  — 

i=:2p-l  i=29-l 

{ [pf  sin  Vi  +  (zf  -  Zj)  cos  Vi]  //i,„(pf ,  pj,  zf  -  Zj)  -  pj  sin  Vi//2,n(pf ,  Pj,  ^f  -  2j)}  , 


(115c) 
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and 


2p+2  2(7+2 

=  E  E  -  z,). 

i=2p-l  j=29-l 


(115(i 


In  the  integrals  (41a-c)  and  (59a-c)  that  define  the  Gk,n  and  Hk,n,  k  —  1, 2, 3  in  (112)  to  115), 
from  (38)  and  (42) 

ll/2 


Rij  = 


(ft  -  Pjf  +  {Zi  -  Zjf  + 


(116) 


and 


{Pi  -  Pjf  +  (4  -  Zjf  +  ^pipj  ) 


11/2 


(117) 


The  two  superscripts  of  the  Z's,  t  or  (f>,  indicate  the  unit  vectors  of  the  testing  functions 
and  expansion  functions,  respectively  (see  (26)  and  (25)).  When  i  =  j,  we  see  from  (116) 
that  Rij  =  0  when  =  0  so  that  the  integrals  defining  the  Gk,n  become  singular.  Hence,  for 
i  =  j,  Rij  is  approximated  by  an  equivalent  distance  [11] 


Rij  = 


(di/4)2  +  4p^sin2(y) 


1/2 


(118) 


Since  Rfj  never  equals  zero  no  special  treatment  using  an  equivalent  distance  is  needed. 

Next,  substituting  the  four-impulse  triangle  function  representation  (96)  of  the  expansion 
functions  in  (69)  followed  by  the  change  in  indices,  we  obtain  the  following  expressions  for 
the  elements  of  the  V  vectors: 

2p+2 

E  T, 

i~2p~\ 

2p+2 

=  E  T7(Cos0‘-(4+l-^J„_,)e>*^•'=-'’•"^  (1195) 

z=2p— 1 
2p+2 

E  +  (119c) 

i=2p-l 


cos  0*"*^  sin  Vi  ( J„+i  —  Jn-i)  +  2y  sin  cos  Uj 


^jkzi  cos  9'^^ 


(119a) 


46 


and 


2p+2 


where 


V4f  =  Jrr+‘  E  r,(^»+i- Jn-Oe'"''"”*'”, 

i=2p— 1 


(119d) 


(119e) 


In  (119a-d)  the  two  superscripts  of  the  Vs  indicate  respectively  the  unit  vectors  of  the 
testing  functions  given  by  (26)  and  the  polarization  of  the  incident  electric  field  given  by 
(11)  and  (12). 

Similarly,  substituting  (96)  in  (85)  followed  by  the  change  in  indices,  we  obtain 

2p+2 


Ti'  [cos 

■1 

Qfar 

sin  Vi 

(«^n+l 

-  Jn-i)  +  2jsin^^“^cosUiJ„] 

,  (120o) 

2p+2 

cos  6^^^ 

=  - 

irf 

E 

Ti'  cos 

+  Jn 

(1206) 

i=2p”l 

2p+2 

>^^jkzi  cos 

^^np 

Txf 

E 

Ti'  sin 

^i('^n+l  d”  Jn—l 

(120c) 

i=2p-l 


and 


2p-f  2 


where 


i=2p~l 


Jn  =  J„(fcpi  sin  0-^“’'). 


(120d) 


(120e) 


In  (120a-d)  the  two  superscripts  of  the  B!s  indicate  the  far  scattered  electric  field  component 
and  the  polarization  of  the  incident  electric  field,  respectively. 

An  alternative  choice  of  the  expansion  and  testing  functions,  closely  related  to  that  of 
Mautz  and  Harrington  which  we  have  used,  is  to  let  the  functions  fi{t)  be  given  by 


h{t)  =  Ti{t) 


(121) 
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instead  of  by  (93),  with  T,(t)  given  by  (92)  as  before.  Then 


j^[pmm  =  i>m(t)+p'(t)r[t)  (122) 

where  r/(t)  is  given  as  before  by  (95).  Again,  using  the  four-impulse  approximation  to 
perform  the  t  and  t'  integrations,  (96)  is  replaced  by 

4 

~  E  Tp+4i-4S{t  -  tp+2i-2)  (123) 

p=l 

and  (98)  is  replaced  by 

j  E  -  V2.-2)  (124) 

with 

"Tp+ii-i  —  fcpp+2t-2Tp+4i— 4,  P  —  1)  ■  ■  ■  )4  (125) 

and 

7j^+4i_4  =  kpp^2i-2'^p4-4i-4  ^  sin  Vp^2i-2Tp4-4i-4,  P  —  Ij  ‘  '  j  4.  (126) 

In  (125)  and  (126)  Tp+4j_4  and  are  given  by  (97)  and  (99)  respectively.  The  expres¬ 

sions  (112)-(120)  for  the  elements  of  the  Z  matrices  and  the  V  and  R  vectors  can  then  be 
used  exactly  as  they  are  simply  by  substituting  T  for  T  and  T'  for  T'.®  It  is  shown  in  Section 
3.6  that  the  alternative  choice  of  the  (121),  must  be  used  instead  of  (93)  if  half-triangle 
functions  are  to  be  used  at  the  beginning  and  end  of  the  discretized  generating  curve  of  the 
BOR. 

®This  efficient  method  of  implementing  the  alternative  choice  of  the  functions  fi{t)  is  due  to  John  Putnam 
[25], 
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2.9  Increasing  the  Efficiency  of  the  Computations 


Referring  back  to  the  general  form  of  the  computations,  (110),  we  observe  that  the  T's 
are  trivial  to  calculate  and  store,  once  and  for  all,  at  the  outset  of  a  calculation  of  the 
far  scattered  field  pattern  for  a  BOR,  whereas  in  contrast  the  F's  are  time  consuming  to 
calculate.  Furthermore  we  observe  that  in  (110)  the  same  F{i,j)  will  contribute  to  the 
calculation  of  several  Qpq.  Hence,  rather  than  using  (110)  to  calculate  Qpq,  it  is  desirable  to 
have  the  limits  of  the  summations  go  from  one  to  the  number  of  segments  in  the  discretized 
generating  curve  of  the  BOR  so  that  there  is  no  repetition  of  the  calculation  of  F{i,j).  The 
basic  question  to  be  answered,  then,  is  given  i  and  j  what  are  the  values  of  p  and  q  for  which 
F{i,j)  contributes  to  Qpq?  An  example  will  help  to  clarify  this  question.  Suppose  i  =  5  and 
j  —  6.  Then  we  can  make  the  following  tabulation 


p  2p  —  I  2p  2p+l  2p  +  2  q  2q  —  \  2q  2q+\  2q  +  2 

1  1  2  3  4  1  1  2  3  4 


10 


10 


and  observe  that  F(5, 6)  contributes  to  Qpq  when  p  is  either  2  or  3  and  q  is  either  2  or  3, 
so  that  F(5, 6)  contributes  to  four  different  Qpq.  In  general  it  is  clear  from  the  limits  of  the 
summations  in  (110)  that  if  i  is  odd  then  p=  (i-  l)/2  and  p=  (i  + 1)/2  whereas  if  i  is  even 
then  p  =  (i  —  2)/2  and  i/2,  and  similarly  for  j  and  q.  Of  course,  when  i  or  j  equals  1  then 
only  one  of  the  two  values  of  p  or  g  is  a  valid  index  for  Qpq  and  similarly  when  i  or  j  is  equal 
to  the  number  of  segments  in  the  discretized  generating  curve. 
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2.10  Choice  of  the  Number  of  Fourier  Modes  for  Expansion  and 
Testing  Functions. 

The  (f)  dependence  of  the  current  given  by  (24)  is  expressed  as  a  summation  from  —N  to  N 
of  the  Fourier  modes  The  value  of  N  can  be  set  equal  to  the  number  of  Fourier  modes 
sufficient  to  represent,  to  the  desired  accuracy,  the  (f)  variation  of  the  tangential  component 
of  the  incident  electric  field  on  the  surface  of  the  BOR.  Let  a  be  the  largest  value  of  p  of  a 
point  (p,  2)  on  the  generating  curve  of  the  BOR.  Then  from  (10),  (11),  and  (12)  it  can  be 
seen  that  the  (j)  variation  of  the  incident  field  along  the  circle  on  the  BOR  corresponding  to 
the  point  (a,  z)  is  given  by 

/((^)  =  or  f{(f))  =  (127) 

For  axial  incidence  0'”*^  =  0  and  /{(f))  equals  cos0  or  sin^  so  that  only  the  modes  are 
needed.  For  oblique  incidence,  we  can  express  /(<?!>)  as  the  Fourier  series 

/(<!>)  =  '£c,e>'^  (128a) 

-N 

with 

C„  =  —  /  (1286) 

0 

from  which,  with  (127)  and  (68a,b), 

I  c„  1=  ^  [Jn+i(fcasin^‘”")  ±  J„_i(/cosin0*"‘=)]  .  (129) 

As  n  increases  beyond  kasin  |  Jn{kasin  |  decreases  rapidly  and  it  suffices  to  choose 

N  =  I  +  M  (130) 


50 


where 


I  =  Int[fcasin^‘""; 


(131) 


and  M  is  the  smallest  integer  for  which 


J^ika  sin^*”*^) 
Ji{ka  sin0'”‘^) 


<  e 


(132) 


with  e  a  small  number  depending  on  the  desired  accuracy.  If  the  value  of  N  given  by  (130) 
is  plotted  as  a  function  of  Intf/casin^®"*^]  the  plot  is  found  to  be  almost  linear.  For  e  =  0.005, 
for  example, 

N  «  Int[l.O4A;asin0'"‘^]  +  7.  (133) 


A  similar  expression 


N  «  Int[A:‘''(asin  6)  +  A] 


(134) 


with  k'^  denoting  a  value  a  few  percent  larger  than  k  was  obtained  by  Yaghjian  [26]  for  the 
reciprocal  problem  of  estimating  the  number  of  angular  modes  needed  to  represent  the  far 
field  of  a  radiator  in  the  9  direction. 


3  RESULTS 

The  analysis  presented  in  Section  2  was  implemented  in  a  FORTRAN  computer  program 
which  was  then  used  to  obtain  numerical  results  for  several  different  BOR’s.  These  results 
will  now  be  described  and  discussed. 


51 


3.1  Spurious  Resonances  for  a  PEC  Sphere 


As  noted  in  the  Introduction,  the  ordinary  EFIE  does  not  have  a  unique  solution  to  the  exte¬ 
rior  scattering  problem  at  frequencies  equal  to  the  resonant  frequencies  of  the  corresponding 
interior  cavity.  Consequently,  numerical  results  obtained  by  solving  the  EFIE  can  have  sig¬ 
nificant  errors  called  spurious  resonances  when  the  electrical  dimensions  of  the  scattering 
object  are  such  as  to  make  the  interior  cavity  resonant.  The  DSEFIE  by  contrast  has  a 
unique  solution  to  the  exterior  scattering  problem  at  all  frequencies.  For  a  sphere,  the  cavity 
problem  has  been  treated  by  Harrington  [27]  who  shows  that  the  condition  for  resonance  is 
ka  =  Unp  for  TE  modes  and  ka  =  for  TM  modes  where  is  the  zero  of  the  Bessel 
function  Jn-\-i/2{^)  and  itjjp  is  the  zero  of 

As  an  example  of  how  the  DSEFIE  eliminates  spurious  resonances  that  may  occur  in 
solving  the  EFIE,  in  Fig.  8  we  show  the  backscatter  RCS  of  a  PEC  sphere  as  a  function  of 
ka  in  the  vicinity  of  the  resonant  value  of  ka  =  6.062,  the  first  zero  of  [x^l'^J^i2]'{x).  It  is 
clearly  seen  that  the  spurious  resonance  displayed  by  the  EFIE  curve  at  6.08  is  eliminated 
by  the  DSEFIE  curve  which  agrees  closely  with  the  exact  Mie  solution  curve. 

3.2  Spurious  Resonance  for  a  PEC  Prolate  Spheroid 

As  a  second  example  of  the  DSEFIE  removal  of  a  spurious  resonance  obtained  using  the 
EFIE,  in  Fig.  9  we  show  the  Eg  bistatic  RCS  pattern  in  the  <f)  =  90°  plane  of  a  prolate 
spheroid  illuminated  by  a  TM  plane  wave.  The  geometry  of  the  spheroid  is  shown  in  Fig. 
10.  The  propagation  vector  of  the  incident  plane  wave  lies  in  the  xz  plane  making  an  angle 
of  0*"'^  =  45°  with  the  z  axis,  and  the  semi-major  and  semi-minor  axes  of  the  spheroid  are 
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Figure  8.  Back  scatter  RCS  of  a  perfectly  conducting  sphere  as  a  function  of  ka  calculated 
by  the  exact  Mie  series,  the  conventional  EFIE,  and  the  DSEFIE.  The  spurious  resonance 
in  the  EFIE  curve  at  ka  =  6.08  is  eliminated  by  the  DSEFIE. 


given  by  ka  =  50  and  kb  =  25,  respectively.  The  large  error  displayed  by  the  EFIE  pattern 
for  6  <  10°  is  eliminated  in  the  pattern  obtained  using  the  DSEFIE  computer  program  and  in 
the  pattern  obtained  with  the  combined-field  BOR  integral  equation  program  CICERO  [21] 
shown  for  reference.  The  discrete  nature  of  the  spurious  resonance  can  be  seen  by  comparing 
Fig.  9  with  Fig.  11  which  shows  the  corresponding  bistatic  RCS  patterns  for  a  spheroid 
with  ka  =  49.8  and  kb  =  24.9.  In  Fig.  11  the  EFIE  pattern  has  no  spurious  resonance. 


3.3  Spurious  Resonances  for  a  PEC  Finite  Cylinder 


Calculations  of  scattering  from  a  finite  cylinder  provide  another  example  of  the  elimination  of 
spurious  resonances  by  the  DSEFIE.  Expressions  for  the  resonant  frequencies  of  a  cylindrical 
cavity  are  derived  in  [28].  In  particular,  the  resonant  frequency  for  the  dominant  TE  mode, 
TEiii,  is  given  by 


Wii 


1.841  / 

1  =  -^=  1  +  2.912— 


^2\  1/2 


y/fioeoR  \ 


(135) 


where  R  and  d  are  the  radius  and  height,  respectively,  of  the  cylinder.  Since  ojy/Ji^  =  k  = 
27r/A,  for  i?  =  d  and  A  =  1  the  radius  and  height  of  the  cylinder  that  supports  the  TEm 
mode  is  R  =  d  =  1.841(3.912)^/^/(27r)  =  0.580.  Figure  12  shows  the  monostatic  E-plane 
RCS  pattern  for  R  =  d  =  .587  as  obtained  by  the  EFIE,  the  CFIE,  and  the  DSEFIE,  with  a 
discretization  of  the  generating  curve  for  the  finite  cylinder  using  20  points/A.  The  striking 
spurious  resonance  displayed  by  the  EFIE  pattern  is  removed  by  both  the  CFIE  and  DSEFIE 
which  are  in  close  agreement  with  each  other.  In  contrast,  the  EFIE,  CFIE,  and  DSEFIE 
monostatic  patterns  for  the  non-resonant  cylinder  R  =  d  =  .52  are  shown  in  Fig.  13.  For 
this  cylinder  there  is  no  spurious  resonance  in  the  EFIE  pattern.  Figure  14  shows  the  EFIE, 
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Figure  1 
at  45°;  k 


CFIE,  and  DSEFIE  curves  for  bistatic  scattering  in  the  E-plane  for  oblique  incidence  at  30° 
on  the  resonant  finite  cylinder,  R  =  d=  .587.  As  in  Fig.  12,  the  spurious  resonance  behavior 
of  the  EFIE  curve  is  removed  by  both  the  CFIE  and  DSEFIE. 

3.4  Bistatic  Scattering  From  the  Cone-Sphere  with  Axial  Illumi¬ 
nation 

The  cone-sphere  is  the  BOR  whose  generating  curve  is  shown  in  Fig.  15.  It  is  determined  by 
two  parameters;  1)  the  semi-vertex  angle  (  of  the  cone,  and  2)  the  radius  a  of  the  spherical 
cap.  The  length  b  of  the  straight-line  segment  of  the  generating  curve  is  given  by  6  =  a/  tan  ^ 
and  the  length  of  the  circular  portion  of  the  generating  curve  is  a(7r/2  -I-  C)-  Also  shown  in 
Fig.  15  is  the  generating  curve  for  the  cone-sphere  dual  surface  separated  by  6  from  the 
original  surface. 

The  use  of  the  DSEFIE  to  calculate  scattering  from  the  cone-sphere  introduces  a  problem, 
not  encountered  in  calculating  scattering  from  the  sphere  or  spheroid,  that  will  be  discussed 
here.  In  Fig.  16  we  have  plotted  bistatic  scattering  in  the  E-plane  from  a  cone-sphere  with 
C  =  5°  and  ka  =  5.9  for  axial  incidence  in  the  — ^  direction,  as  calculated  by  the  CFIE  and 
DSEFIE.  The  separation  S  of  the  dual  surface  from  the  original  surface  is  5  =  A/4.  It  is  seen 
that  there  are  small  oscillations  in  the  central  portion  of  the  DSEFIE  curve  that  are  not 
present  in  the  combined-field  integral  equation  pattern  obtained  using  the  computer  code 
CICERO.  To  investigate  the  source  of  these  oscillations,  in  Fig.  17  we  have  plotted  the  CFIE 
RCS  pattern  of  Fig.  16  along  with  the  two  “components”  of  the  DSEFIE  pattern:  the  EFIE 
pattern  obtained  with  the  original  cone-sphere  surface  alone,  and  the  DSEFIE  pattern 
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Figure  16.  E-plane  bistatic  RCS  patterns  for  axial  incidence  in  the  —z  direction  on  the 

cone-sphere,  C  =  5°,  A:a  =  5.9;  5  =  A/4  for  the  DSEFIE  curve.  Note  the  small  oscillations  in 
the  central  portion  of  the  DSEFIE  curve. 


obtained  with  the  dual  surface  aloneJ  It  is  seen  that  while  the  RCS  curve  obtained  with  the 
EFIE  and  the  original  cone-sphere  surface  is  very  close  to  that  obtained  with  the  CFIE,  the 
RCS  curve  obtained  using  the  DSEFIE  with  the  boundary  condition  enforced  on  the  dual 
surface  alone  differs  greatly  from  the  other  two  curves.  With  the  dual-surface-only  pattern 
so  much  in  error,  the  DSEFIE  cannot  be  expected  to  give  accurate  results,  since  the  DSEFIE 
calculations  will  be  very  sensitive  to  the  value  of  the  combination  parameter,  a.  For  further 
understanding,  the  currents  on  the  cone-sphere  surface  were  printed  out  and  it  was  found 
that  when  the  zero-tangential  electric  field  boundary  condition  was  enforced  on  the  dual 
surface  alone,  the  currents  close  to  the  vertex  of  the  cone  were  several  orders  of  magnitude 
larger  than  the  currents  close  to  the  vertex  when  the  zero-tangential  electric  field  boundary 
condition  was  enforced  on  the  original  cone-sphere  surface  alone.  The  discrepancy  between 
the  DSEFIE  and  CFIE  RCS  patterns  in  Fig.  16  is  therefore  due  to  a  problem  in  calculating 
scattering  from  the  cone-sphere  when  a  zero-tangential  electric  field  is  enforced  on  the  dual 
surface  alone. 

A  partial  explanation  for  why  this  problem  occurs  appears  to  be  the  following.  Unlike 
for  the  sphere  or  spheroid  where  corresponding  points  on  the  original  and  dual  surfaces  are  a 
specified  distance  5  apart,  for  the  cone-sphere  points  close  to  the  vertex  can  be  considerably 

further  apart  than  5  from  their  corresponding  points  on  the  dual  cone-sphere  surface.  For 

^When  we  refer  to  the  “components”  of  the  DSEFIE  pattern,  we  do  not  mean  to  imply  that  the  total 

DSEFIE  pattern  is  a  linear  combination  of  the  EFIE  solution  and  the  DSEFIE  solution  obtained  with  the 
zero-tangential  electric  field  boundary  condition  enforced  on  the  dual  surface  alone.  This  is  not  the  case. 
Loosely  speaking,  however,  there  is  no  harm  in  referring  to  these  two  limiting  forms  of  the  DSEFIE  solution 
as  “components” ,  especially  as  considerable  insight  into  the  DSEFIE  solution  can  be  obtained  by  so  doing. 
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example,  the  vertex  point  on  the  dual  cone-sphere  surface  is  a  distance  5/sin(^  from  the 
vertex  of  the  original  cone-sphere.  For  C  =  5°  and  5  =  A/4  this  distance  is  almost  3A.  The 
boundary  condition  of  zero-tangential  electric  field  on  the  dual  surface  is  then  apparently 
insufficient  to  guarantee  that  the  currents  close  to  the  tip  of  the  cone-sphere  will  be  correct. 
It  should  be  noted,  however,  that  this  problem  is  specifically  associated  with  the  tip  of  the 
cone  and  is  not  simply  a  result  of  corresponding  points  on  the  original  and  dual  surface 
varying  in  distance  from  one  another.  For  the  finite  cylinder  generating  curves  shown  in 
Fig.  4,  for  example,  the  distance  between  the  corresponding  points  on  the  corners  of  the 
original  and  dual  surface  generators  is  greater  by  a  factor  of  \/2  than  the  distance  between 
the  corresponding  points  on  the  midpoints  of  the  horizontal  or  vertical  line  segments,  but 
there  is  no  problem  in  calculating  the  RCS  pattern  when  the  zero-tangential  electric  field 
boundary  condition  is  enforced  on  the  dual  surface  alone. 

The  problem  we  have  described  cannot  be  solved  simply  by  placing  the  dual  surface 
considerably  closer  than  A/4  to  the  original  surface.  Even  a  separation  of  A/64  does  not 
satisfactorily  resolve  the  difficulty.  What  can  be  done,  however,  in  addition  to  placing  the 
dual  surface  closer  to  the  original  surface  is  to  insert  an  auxiliary  dual  surface  close  to  the 
vertex  of  the  cone-sphere  as  shown  in  Fig.  18  in  which  the  line  segment  that  generates  this 
auxiliary  dual  surface  is  obtained  by  drawing  a  line  from  the  origin  to,  say,  the  midpoint  of 
the  line  from  the  vertex  of  the  dual  cone-sphere  normal  to  the  original  cone-sphere  generator. 
The  point  on  this  generating  line  segment  of  the  auxiliary  dual  surface  that  corresponds  to 
the  vertex  of  the  original  cone-sphere  is  taken  a  small  distance,  say  A/64,  along  this  line.  In 
Fig.  19  we  show  the  bistatic  RCS  patterns  for  the  same  cone-sphere  as  in  Fig.  16  obtained 
using  the  DSEFIE  with  an  enhanced  dual  surface  with  S  =  A/32  for  the  main  dual  surface. 
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The  small  oscillations  present  in  the  DSEFIE  pattern  of  Fig.  16  are  completely  removed. 
Also  shown  in  Fig.  19  is  the  pattern  obtained  with  the  DSEFIE  when  only  the  enhanced  dual 
surface  is  employed.  In  striking  contrast  to  the  DSEFIE  pattern  of  Fig.  17  the  dual-surface- 
only  pattern  of  Fig.  19  is  now  very  close  to  the  CFIE  pattern.  To  demonstrate  that  the 
great  improvement  obtained  with  the  enhanced  dual  surface  is  the  result  of  both  employing 
the  auxiliary  dual  surface  in  the  vicinity  of  the  cone-sphere  vertex  and  moving  the  main  dual 
surface  closer  to  the  original  surface,  in  Fig.  20  we  have  plotted  two  additional  patterns 
along  with  the  CFIE  pattern.  One  pattern  corresponds  to  employing  the  same  auxiliary 
dual  surface  but  with  the  main  dual  surface  separated  by  A/4  from  the  original  surface,  and 
the  other  pattern  corresponds  to  simply  moving  the  main  dual  surface  close  to  the  original 
cone-sphere  surface  by  taking  S  =  A/32  but  not  employing  the  auxiliary  dual  surface.  It  is 
seen  that  employing  the  auxiliary  dual  surface  while  keeping  the  main  dual  surface  separated 
by  A/4  from  the  original  surface  results  in  considerable  improvement  compared  to  the  dual¬ 
surface-only  pattern  of  Fig.  17  though  not  as  much  as  that  obtained  when  the  main  dual 
surface  is  moved  closer  and  the  auxiliary  dual  surface  employed.  Merely  moving  the  main 
dual  surface  close  to  the  original  surface  without  employing  the  auxiliary  dual  surface  still 
results  in  a  very  poor  RCS  pattern,  however.  For  larger  semi- vertex  angles  C  the  problem 
we  have  described  in  this  section  is  less  severe  but  is  still  present.  For  C  =  45°,  for  example, 
a  separation  of  (5  =  A/32  suffices  to  resolve  the  problem  without  the  necessity  of  inserting 
an  auxiliary  dual  surface  close  to  the  cone  tip. 
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3.5  Cone-Sphere  Monostatic  Scattering 


In  this  section  we  discuss  scattering  for  oblique  incidence  on  the  cone-sphere  in  contrast 
to  the  previous  section  where  axial  incidence  was  considered.  For  axial  incidence,  only  the 
first  order  Fourier  mode  is  required  to  describe  azimuthal  variation.  For  oblique  incidence, 
however,  the  zeroth  order  Fourier  mode  as  well  as  higher  order  modes  are  required  in  addition 
to  the  first  order  mode.  It  will  be  seen  that  the  presence  of  the  zeroth  order  rriode  in  particular 
introduces  complications  not  encountered  for  axial  incidence. 

Figure  21  shows  the  monostatic  pattern  in  the  E-plane  for  the  narrow  tip-angle  cone- 

sphere,  C  =  5°,ka  =  5.9,  treated  in  Section  3.4  as  calculated  by  the  CFIE®,  EFIE,  and 

DSEFIE,  using  Fourier  modes  0  through  6.  The  DSEFIE  calculation  is  performed  with  the 

enhanced  dual  surface  described  in  Section  3.4,  with  a  separation  of  <5  =  A/32  between  the 

original  cone-sphere  surace  and  the  main  dual  surface.  Figure  22  shows  the  corresponding 

curves  in  the  H-plane.  It  is  seen  that  there  is  a  significant  difference  between  the  values  of 

the  E-plane  RCS’s  for  6  in  the  vicinity  of  15°,  while  the  H-plane  patterns  exhibit  only  small 

differences  in  a  low  dB  region  of  the  patterns.  In  Fig.  23  we  show  the  E-plane  monostatic 

patterns  when  modes  1  through  6  are  used,  and  observe  that  there  is  virtually  no  difference 

between  the  CFIE,  EFIE,  and  DSEFIE  patterns.  In  contrast,  the  CFIE,  EFIE,  and  DSEFIE 

patterns  for  the  zeroth  order  mode  only,  shown  in  Fig.  24,  display  significant  differences  in 

the  vicinity  of  6  =  15°.  The  results  shown  so  far  thus  demonstrate  that  the  diflferences  in  the 

E-plane  monostatic  patterns  plotted  in  Fig.  21  are  attributable  to  differences  in  the  zeroth 

order  mode  patterns.  Furthermore,  since  the  E-plane  and  H-plane  zeroth  order  mode  pat- 
*Unless  otherwise  noted,  when  we  refer  to  the  CFIE  it  is  assumed  that  the  combination  parameter  oq 

(see  (3))  is  set  equal  to  0.5,  thereby  giving  equal  weight  to  the  EFIE  and  MFIE  components  of  the  CFIE. 
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Figure  22.  H-plane  monostatic  RCS  patterns  for  the  cone-sphere,  C  =  5°,  ka  —  5.9;  Fourier 
modes  0  through  6.  Enhanced  dual  surface  used  for  the  DSEFIE  calculation  with  S  =  A/32. 
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terms  are  due  solely  to  the  t-directed  and  (/(-directed  currents,  respectively,  it  is  the  difference 
in  accuracy  with  which  the  fields  radiated  by  the  zeroth  order  mode  t-directed  currents  are 
calculated  by  the  CFIE,  EFIE,  and  DSEFIE,  that  accounts  for  the  differences  we  have 
observed. 

To  establish  which  one  of  the  CFIE,  EFIE,  and  DSEFIE  patterns  of  Fig.  24  can  be 
regarded  as  correct,  we  begin  by  considering  the  discrepancy  between  the  CFIE  and  EFIE 
patterns.  Since  it  is  known  that  the  MFIE  is  ill  suited  to  calculations  involving  scatterers 
with  surfaces  separated  by  much  less  than  a  wavelength,  as  is  the  case  with  narrow  tip-angle 
cones,  it  can  be  conjectured  that  the  discrepancy  between  the  CFIE  and  EFIE  patterns  is 
due  to  the  MFIE  “component”  that  enters  into  the  CFIE  in  combination  with  the  EFIE, 
thereby  degrading  the  accuracy  of  the  EFIE  solution.®  In  Fig.  25  we  have  plotted  the  EFIE 
and  MFIE  zeroth  order  mode  El-plane  monostatic  patterns.  We  observe  that  the  MFIE  is 
considerably  at  variance  with  the  EFIE  not  only  in  the  vicinity  of  0  =  15°  but  for  almost 
the  entire  range  of  6  from  0°  to  90°.  To  verify  that  it  is  indeed  the  MFIE  pattern  that  is 
incorrect,  in  Fig.  26  we  have  plotted  the  EFIE  pattern  of  Fig.  23  along  with  the  MFIE 
pattern  as  recalculated  with  a  density  of  400  points/ A  for  the  first  wavelength  from  the  tip 
of  the  cone  where  the  t-currents  vary  most  rapidly,  and  a  density  of  40  points/ A  elsewhere. 
The  MFIE  pattern  is  now  closer  to  the  EFIE  pattern,  though  it  still  displays  a  discrepancy 
of  almost  2  dB  in  the  near  vicinity  of  0  =  15°.  Also  shown  in  Fig.  26  is  the  EFIE  pattern 

®When  we  refer  to  the  MFIE  and  EFIE  “components”  of  the  CFIE,  we  do  not  mean  to  imply  that  the 
CFIE  solution  is  a  linear  combination  of  the  MFIE  and  EFIE  solutions.  There  is  no  harm,  however,  in 
thinking  of  these  two  limiting  forms  of  the  the  CFIE  solution  as  “components”  since  considerable  insight 
into  the  CFIE  solution  can  be  obtained  by  so  doing. 
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Figure  25.  E-plane  monostatic  EFIE  and  MFIE  RCS  patterns  for  the  cone-sphere, 
ka  =  5.9;  Fourier  mode  0  only. 
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as  calculated  with  a  400  point/ A  density  close  to  the  tip.  It  is  seen  that  there  is  almost  no 
change  in  the  EFIE  pattern.  We  therefore  feel  justified  in  regarding  the  EFIE  pattern  as 
correct. 

Now  the  CFIE  was  developed  to  remove  spurious  resonances  resulting  from  non-uniqueness 
of  the  MFIE  and  EFIE  at  internal  cavity  resonant  frequencies.  For  non-resonant  scatterers, 
such  as  the  particular  cone-sphere  we  are  concerned  with,  the  CFIE  gives  correct  results 
provided  that  the  MFIE  and  EFIE  solutions,  separately,  are  correct.  Here,  however,  as  we 
have  seen,  the  MFIE  solution  of  Fig.  25  that  enters  into  the  CFIE  solution  shown  in  Fig. 
24  is  highly  inaccurate.  Combining  an  inaccurate  MFIE  solution  with  an  accurate  EFIE 
solution  can  only  degrade  the  accuracy  of  the  EFIE  solution,  as  witness  the  discrepancy 
between  the  CFIE  and  EFIE  patterns  of  Fig.  24.  If  the  CFIE  pattern  is  recalculated  using 
a  density  of  400  points/ A  close  to  the  cone  tip,  then  the  CFIE  pattern  is  much  closer  to  the 
EFIE  pattern,  as  shown  in  Fig.  27.  This  is  as  expected  since,  as  we  have  seen,  the  accuracy 
of  the  MFIE  component  of  the  CFIE  is  thereby  considerably  increased.  The  problem  we 
have  encountered  here  with  the  CFIE  corroborates  the  finding  of  Wood  and  Hill  [23]  that 
the  value  of  the  combination  parameter  used  in  the  CFIE  can  strongly  influence  the  rate 
of  convergence  of  the  solution  as  the  density  of  points/ A  is  increased.  An  accurate  EFIE 
(ao  =  1-0)  solution  can  be  obtained  here  with  a  40  points/ A  density,  but  a  value  of  ao  =  0.5 
(equal  weighting  of  the  EFIE  and  MFIE)  requires  a  density  of  approximately  400  points/A 
to  yield  acceptable  results. 

Having  established  the  accuracy  of  the  EFIE  zeroth  order  mode  E-plane  monostatic 
pattern,  how  is  the  discrepancy  shown  in  Fig.  24  between  the  DSEFIE  and  EFIE  solutions 
to  be  explained?  As  with  the  CFIE,  employing  a  density  of  400  points/ A  for  the  first  wave- 
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length  from  the  cone  vertex  and  a  density  of  40  points/ A  elsewhere  gives  a  DSEFIE  pattern 
very  close  to  the  EFIE  pattern  as  shown  in  Fig.  28.  It  is  difficult  to  be  sure  without 
extensive  further  investigation  why  this  high  density  of  points  close  to  the  cone  tip  is  needed 
to  obtain  a  highly  accurate  solution  with  the  DSEFIE.  Since  the  EFIE  by  itself  gives  an 
acceptable  pattern  for  a  density  of  40  points/A,  it  is  clearly  the  calculation  of  the  dual 
surface  component  of  the  DSEFIE  that  requires  a  high  density  of  points.  Furthermore,  since 
a  density  of  80  points/A  gave  acceptable  DSEFIE  results  for  the  axial  incidence  bistatic 
scattering  discussed  in  Section  3.4,  it  appears  that  the  rapid  variation  of  the  zeroth  order 
mode  t-directed  currents  in  combination  with  some  aspect  of  the  dual  surface  integration 
procedure,  possibly  the  treatment  of  the  self-term,  is  responsible  for  the  difficulty  with  the 
dual  surface  calculation.  The  essential  thing  to  stress,  however,  as  we  did  in  connection  with 
the  CFIE,  is  the  importance  of  numerical  experimentation  before  accepting  the  validity  of 
surface  integral  scattering  calculations,  especially  if  high  accuracy  is  desired. 

Since,  as  stated  in  the  Introduction,  the  rationale  for  developing  the  DSEFIE  has  been  the 
importance  of  having  two  independent  BOR  surface  integral  equation  scattering  formulations 
and  computer  codes,  and  since,  as  we  have  also  noted  in  the  Introduction,  a  DSMFIE 
formulation  already  exists,  it  is  of  interest  to  see  how  well  the  DSMFIE  calculates  the  zeroth 
order  mode  monostatic  E-plane  pattern  for  the  narrow  tip-angle  cone  considered  here.  In 
Fig.  29  we  show  the  DSMFIE  pattern  calculated  with  a  density  of  400  points/A  close  to 
the  cone  tip  and  40  points/A  elsewhere,  along  with  the  EFIE  pattern  of  Fig.  24,  and  the 
DSEFIE  pattern  of  Fig.  28.  A  separation  of  5  =  A/4  between  the  dual  surface  and  the 
original  cone-sphere  surface  is  used  for  the  DSMFIE  calculation  (the  auxiliary  dual  surface 
used  for  the  DSEFIE  cone-sphere  calculations  is  not  useful  for  the  DSMFIE  calculations). 
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Figure  28.  E-plane  monostatic  EFIE  and  DSEFIE  RCS  patterns  for  the  cone-sphere,  C  = 
ka  =  5.9;  Fourier  mode  0  onl3\  The  DSEFIE  pattern  is  obtained  using  a  400  point /A  erid 
n  the  t  direction  for  a  wavelength  from  the  cone-tip  and  a  density  of  40  points/A  elsewhere 


The  DSMFIE  pattern  is  in  error  by  1.6  dB  at  0  =  13°  and  there  are  additional  discrepancies 
in  the  location  of  peaks  and  nulls.  We  conclude  that  for  calculations  involving  scatterers 
with  narrow  tips,  the  DSEFIE  is  significantly  more  reliable  than  the  DSMFIE  (see,  however. 
Footnote  3). 


3.6  Half- Triangle  Functions 


In  Section  2.8  the  expansion  and  testing  functions  were  defined  in  terms  of  triangle  functions 
Ti{t)  shown  in  Fig.  6.  The  Ti{t)  overlap  as  shown  in  Fig.  7  so  that,  for  example,  =  1 

at  ^21-1  where  Ti{t)  =  0,  and  =  1  at  t^i+i  where  Ti{t)  =  0,  etc.  At  the  lower  and  upper 

poles  of  the  discretized  generating  curve  there  is  no  overlap,  however,  since  the  first  triangle 
function  is  Ti(t)  which  is  0  at  the  lower  pole,  and  the  last  triangle  function  is  r(p_  3)/2(0 
which  is  0  at  the  upper  pole.  In  certain  applications  it  may  be  desirable  to  have  overlapping 
at  the  lower  and  upper  poles  of  the  generating  curve.  This  can  be  achieved  by  employing 
the  upper  half  of  the  triangle  function  To{t)  at  the  lower  pole  where,  from  (92), 


ro{t)  =  { 


0, 

and  adding  the  lower  half  of  the  triangle  function  r(p_3)/2+i(t)  at  the  upper  pole  where 


(136) 


T{P-Z)/2+li't)  -  S 


0, 


t  <  tp_2,  t  >  t*p 


t  ^P— 2  f*  <'#<'#* 

I  dp_2  +  dp-i’ 


(137) 


If  half-triangle  functions  are  employed  at  the  beginning  and  end  of  the  generating  curve, 
however,  it  is  important  to  define  the  expansion  and  testing  functions  fi{t)  by  (121)  rather 
than  by  (93)  if  p  0  at  the  lower  or  upper  pole  of  the  BOR  generating  curve.  If,  for  example. 
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(93)  is  used  to  define  fo{t)  in  conjunction  with  (136)  and  p  — >  0  as  t  0,  fo{t)  becomes 
singular  at  p  =  0.  Even  if  the  actual  value  of  t  =  0  is  not  used  because  of  the  four-impulse 
approximation  (96),  calculations  based  on  (93)  and  employing  half-triangle  functions  can 
yield  erroneous  results. 

As  an  example  of  the  problems  that  can  result  by  employing  half-triangle  functions  along 
with  the  expansion  and  testing  functions  of  (93),  in  Fig.  30  we  show  the  E-plane  RCS 
pattern  of  a  sphere  with  ka  =  5.5  obliquely  illuminated  by  a  TM  plane  wave  at  an  angle  of 
incidence  of  45°.  Along  with  the  Mie  series  pattern  shown  for  reference  we  show  the  DSEFIE 
patterns  obtained  with  the  basis  functions  of  (93)  with  and  without  half-triangle  functions, 
and  the  pattern  obtained  with  the  basis  functions  of  (121)  using  half-triangle  functions.  It  is 
seen  that  the  pattern  obtained  using  half-triangle  functions  and  the  basis  functions  of  (93) 
differs  greatly  from  the  other  three  patterns.  This  discrepancy  does  not  occur  with  the  H- 
plane.  Further  calculations  demonstrate  that  the  discrepancy  that  occurs  using  half-triangle 
functions  along  with  the  basis  functions  of  (93)  is  due  to  the  zeroth  order  Fourier  mode  as 
shown  in  Fig.  31.  No  discrepancy  occurs  with  the  higher-order  Fourier  modes.  Although 
not  an  explanation  it  is  noted  that  for  the  zeroth  order  mode  the  E-plane  pattern  is  due 
solely  to  the  t-directed  currents  and  not  to  the  (^-directed  currents  while  the  H-plane  patten 
is  due  solely  to  the  (^-directed  currents  on  the  sphere. 

4  SUMMARY 

The  purpose  of  this  report  is  to  provide  for  the  first  time  a  detailed  analysis  and  solution  of 
the  problem  of  determining  scattering  from  a  PEC  BOR  using  the  DSEFIE,  and  to  show 
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results  obtained  with  a  computer  code  based  on  this  analysis.  The  DSEFIE  is  an  extension 
of  the  EFIE  in  which  the  boundary  condition  that  the  tangential  electric  field  vanish  on  the 
surface  of  a  PEC  object  is  satisfied  not  only  on  the  actual  surface  but  also  on  a  fictitious 
surface  lying  inside  the  actual  surface  and  separated  by  a  small  distance  from  it.  Scatter¬ 
ing  solutions  obtained  using  the  DSEFIE  are  free  from  the  spurious  resonances  that  can 
seriously  degrade  the  accuracy  of  solutions  obtained  using  the  conventional  EFIE  or  MFIE 
for  3-D  scatterers  whose  linear  dimensions  exceed  a  wavelength.  These  spurious  resonances 
are  a  consequence  of  the  fact  that  the  EFIE  and  MFIE  fail  to  produce  a  unique  solution 
for  the  current  on  a  PEC  scatterer  at  frequencies  equal  to  the  resonant  frequencies  of  the 
interior  cavity  formed  by  the  surface  of  the  scatterer.  Although  the  widely  used  CFIE  also 
eliminates  spurious  resonances,  it  is  important  to  have  available  two  independent  surface 
integral  formulations  that  can  be  applied  and  solved  efficiently.  Another  surface  integral 
equation  solution  free  from  spurious  resonance  that  has  been  successfully  programmed  and 
tested,  the  DSMFIE,  is  also  independent  of  the  CFIE,  but  it  does  not  yield  accurate  results 
for  conductors  containing  narrow-angle  wedges  or  tips. 

In  this  report  the  DSEFIE  is  formulated  in  terms  of  the  vector  and  scalar  potentials,  A 
and  <I>,  and  then  solved  by  the  Galerkin  form  of  the  method  of  moments.  It  is  noteworthy 
that  the  device  of  Mautz  and  Harrington  used  to  transfer  the  differential  operator  on  $ 
to  the  testing  functions  and  so  maintain  the  singularity  of  the  kernel  equal  to  that  of  the 
free-space  Green’s  function,  cannot  be  used  when  the  observation  point  lies  on  the  dual 
surface  rather  than  on  the  original  surface.  As  a  result  the  elements  of  the  Z  matrices  for 
the  DSEFIE  have  a  more  singular  kernel  than  that  of  the  free-space  Green’s  function.  This, 
however,  does  not  appear  to  create  numerical  difficulties  since  the  distance  between  points 
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on  the  original  surface  and  points  on  the  dual  surface  never  goes  to  zero.  The  basis  and 
testing  functions  are  given  as  products  of  the  Fourier  modes  used  to  express  the  azimuthal 
dependence  of  the  currents,  and  triangle  functions  are  used  to  express  the  dependence  of  the 
currents  on  the  BOR  generating  curve  parameter  t.  A  four-impulse  approximation  to  the 
triangle  functions  is  used  which  enables  the  integrations  with  respect  to  t  to  be  performed 
in  closed  form.  Detailed  expressions  for  the  elements  of  the  Z  matrices  and  the  V  column 
vectors  that  appear  in  the  Galerkin  matrix  formulation  of  the  DSEFIE  are  then  given  as 
summations  that  can  be  readily  evaluated  by  computer.  The  efficiency  of  evaluating  these 
expressions  can  be  greatly  increased  by  using  a  change  of  indexing  described  in  the  report. 
Calculations  performed  with  a  computer  program  of  the  DSEFIE  solution  demonstrate  the 
removal  of  spurious  resonances  that  appear  in  calculations  of  the  RCS’s  of  spheres,  spheroids, 
and  finite  cylinders  made  with  the  conventional  EFIE.  Cone-sphere  calculations  show  that 
application  of  the  DSEFIE  to  objects  containing  narrow  tips  requires  careful  placement  of 
the  dual  surface  in  the  vicinity  of  the  tips  and  sometimes  an  increased  density  of  grid  points 
in  the  t  direction  a  wavelength  or  so  from  the  tip.  When  these  precautions  are  taken,  the 
DSEFIE  yields  significantly  more  accurate  solutions  than  the  DSMFIE  for  scatterers  with 
narrow  tips. 
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